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A grid of NLTE corrections for magnesium and calcium in 
late-type giant and supergiant stars: application to Gaia 
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ABSTRACT 

We investigate NLTE effects for magnesium and calcium in the atmospheres of 
late-type giant and supergiant stars. The aim of this paper is to provide a grid 
of NLTE/LTE equivalent width ratios W/W* of Mg and Ca lines for the follow- 
ing range of stellar parameters: T^s S [3500, 5250] K, log g G [0.5, 2.0] dex and 
[Fe/H] € [—4.0, 0.5] dex. We use realistic model atoms with the best physics avail- 
able and taking into account the fine structure. The Mg and Ca lines of interest are in 
optical and near IR ranges. A special interest concerns the lines in the Gaia spectro- 
graph (RVS) wavelength domain [8470, 8740] A. The NLTE corrections are provided 
as function of stellar parameters in an electronic table as well as in a polynomial form 
for the Gaia/KVS lines. 

Key words: line: formation, radiative transfer, stars: late-type, stars: atmospheres, 
stars: abundances. 



1 INTRODUCTION 

Individual stars of different ages carry information that al- 
low astronomers to detail the sequence of events which built 
up the Galaxy. Calcium (Ca) and magnesium (Mg) are key 
elements for the understanding of these events since these 
a-elements are mainly produced by supernovae of type II. 
Therefore, they are both important tracers to derive trends 
of chemical evolution of the Galaxy. These two elements are 
easily observed in spectra of late-type stars through strong 
or weak lines. Among them, there are the Mgi b green 
triplet, the Cai red triplet and the Can IR triplet (CaT) 
lines which are measurable in the most extreme metal-poor 
stars IjChristlieb et al.ll2002l ) even at moderate spectral res- 
olutions. 

The lines of these a-elements provide diagnostic tools 
to probe a wide range of astrophysical mechanisms in stellar 
populations. The study of a-element synthe sis in halo metal- 
poor stars constrains the rate of SNe II (jNakamura et al.l 
11993 ). The Mg element is a reliable reference of the early 
evolutionary time scale of the Galaxy (e.g. Cayrcl ct al. 
l2004lGehren et akllJOO^ and lAndrievskv et al.lboiO. '). More- 
over, the Cai triplet at AA 6102, 6122 and 6162, can be 
used as a gravity in dicator in detailed analyses of dwarf 
and subgiant stars (Cayrel et al. 1996). Two Ca lines of 
different ionization stages can also be used to constraint 
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surface gravity in extremely met al-poor stars (Cai 4226 A 
and Can 8498 A) as shown by iMashonkina et all (|20o3). 
Moreover, the line strengths of CaT lines can be used as 
metallicit y indic ator for stellar population in distant galax- 
ies (jldiart et al.i il997i ) or to search for the most extreme 
metal-po or stars in the vicinity o f nearby dwarf spheroidal 
galaxies (jStarkenburg et al.ll201CI ). The CaT lines are also 
used to probe stellar activity by the measure of the cen- 
tral line depressio n as done by lAndretta et al.l ((20051) and 
iBusa et all (|2007l l. 

It is well known that the assumption of Local Ther- 
modynamic Equilibrium (LTE) is the main source of un- 
certainties (after the quality of the atomic data). In his re- 
view, Asplund (2005) summarized the need to take Non- 
LTE (NLTE) effects into account to get reliable and precise 
chemical abundances. Such effects are crucial in the atmo- 
spheres of metal-poor stars due to several mechanisms: the 
overionization in the presence of a weak opacity, the pho- 
ton pumping and the photon suction. However, since the 
NLTE computations are still time consuming as soon as 
we deal with realistic model atoms, most of the studies are 
still done under the LT E assumption for large sur veys (e.g.. 



Edvardsson et al.l 19931 for lat e-type dwarf stars, iLiu et abl 



2OI0I and iDe Silva et al.l I2OIII for late-type giant and sub- 



giant stars). 

Several studies of NLTE effects on Mgi lines were per- 
form ed for stars of v arious metalli cit ies IIMashonkina et al.l 
1 19961 : IZhao et all Il998l , |200GI : ICratton et all Il999l : 
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Shimanskava et al.1 [2OO0I: lldiart fc ThevenirJ I2OO0I : 
Mishenina et all I2OO4I : lAndrievskv et al.1 I2OI0I '). All stud- 



ies showed an increase of NLTE effects with decreasing 
metallicity which imply corrections up to +0.2 dex. The 
first attempt t o build a reali st ic mo del atom of Mgi was 
performed by ICarlsson et al.1 () 19921 ') to explain the Mgi 
solar emission lines near to 12 /im. They demonstrated 
the photospheric nature of these lines by using the photon 
suction mechanism which required a model atom with 
Rydberg states to ensure the coUisional-dominated coupling 
with the ionization stage. The mo s t com plete Mg model 
atom was built by IPrzvbilla et al.l (|200ll ) with 88 energy 
levels for Mgi and 37 levels for Mgii. However, they did not 
take into account the fine structure of the triplet system. 
The first Mg i model atom that i nclude s the fine structure 
was built by lldiart fc TheveninI (|200Clf ) who applied it to 
the determination of the Mg abundance for more than 
two hundred late-type star s. We can mention a n update 
of the Mgi model ato m of ICarlsson et al.l (|l992l ) done by 
ISundgvist et al.l (12003) in order to reproduce IR emission 
lines in K giants, taking into account the electronic collisions 
between states with the principal quantum number n ^ 15. 

NLTE eff ects on Ca were investigated by 

|j0rgensen et al.l l|l992f l who gave a relation between 
equivalent widths of CaT lines as functions of the surface 
gravities and the ef fectiv e temperatures. As for the Mgi, 
lldiart fc TheveninI (|2000|) included the fine structure in 
their model atom of Cai. They combin ed these results to 
the NLTE Fe abundance predictions (|Thevenin fc Idiarj 
[1999I) and performed an NLTE analysis of the same set 
of late-type dwarf and subgiant stars. They concluded 
that the [Ca/Fe] remained mostly unchanged but some 
parallel structures appeared in the NLTE [Ca/Fe] vs 
[Fe/H] diagram. The mo st complete Ca atom was done by 
iMashonkina et al.l (|2007|) who built a model atom through 
two ionization stages with 63 energy levels for Cai and 37 
levels for Can, including fine structure only for the first 
states. They used the most recent and acc urate atomic data 
from the Opacity Project ljSeatonlll994l 'l. They concluded 
that NLTE abundance corrections can be as large as 0.5 dex 
(for the Cai 4226 A resonance line) in extremely met a l-poor 
stars. Using these NLTE corrections, lAsplund et al.l (|2009|) 
found good agreement of solar Cai and Can abundances. 



In the view of the ESA Gaia mission (jWilkinson et al.l 

[2OO5I), it is important to investigate such NLTE abundance 
corrections in order to distinguish different stellar popula- 
tions. The onboard spectrograph, nam ed RVS (Radial Ve- 
locity Spectrometer, iKatz et al.l |200J) which was designed 
primarily to get the radial velocities can also be used in 
principle to extract individual chemical abundances. The 
spectral window considered was selected in the near-infrared 
[8470-8740] A because of the presence of the CaT lines 
which will be visible even in the most metal-poor stars and 
at a considered low resolution [R ~ 11500 at best). The 
counterpart ground-ba sed observations is the RAVE survey 
ISteinmetz et al.ll2006l ) which will provide about 1 million 
of spectra. These two surveys will analyze a representative 
sample of stars throughout the Galaxy, allowing an unprece- 
dented kinematic and chemical study of the stellar popula- 
tions. The huge amount of stars will prevent us to analyze 
them one by one. It is necessary to have precalculated correc- 



tion tables that can be used for automatic stellar parameter 
determination and this is the aim of this work. 

The structure of the paper is as follows: in Sect. 2, we 
present the model atoms of Mgi, Cai and Can. In Sect. 3, 
we test these models by comparing NLTE synthetic lines 
with several solar (Sect. 3.1) and Arcturus (Sect. 3.2) lines 
of Mgi and Can in the Gaia/RVS region. Then, in Sect. 4.1, 
we present a grid of NLTE corrections of equivalent widths 
of spectral lines, selected in Sect. 4.2, useful for the Gaia 
mission and the RAVE survey. The behaviours of Mgi, Cai 
and Can NLTE lines, the comparisons with previous works, 
the estimate of uncertainties on atomic data, and the poly- 
nomial fits of the Gaia/RVS lines are discussed in Sect. 4. 



2 Mgi, Cai AND Can MODEL ATOMS 

The construction of a model atom consists in collecting and 
computing atomic data which are: 

• the excited states of the atom characterized by their 
energy levels and their statistical weights; 

• the radiative transitions characterized by the oscillator 
strength values / for the bound-bound transitions, and by 
the photoionization cross-sections for the bound-free ones; 

• the collisional transitions with electrons and neutral hy- 
drogen characterized by the collisional strength Q.; 

• the elastic radiations and collisions with neutral hydro- 
gen characterized by the radiative (7rad) and the collisional 
(7coii) damping factors. 



2.1 Construction of model atoms 

Energy levels. The e ne rgy levels come from the 

Ni siB (iR alchenk o et al] |2008| ) and the TopBas((!| 

line atomic databases. 



(|Cunto fc M cndoza Il992l ) 



on 



When theoretical levels are present in the TopBase but 
absent in the nist, we add them in the model atoms by 
shifting them with respect to the nist ionization energies. 
We consider the fine structure of energy levels as long as 
they exist in nist database. If some fine levels are missing 
for the highest levels in atomic databases, we compute them 
as specified in the next sub-sections. 

Radiative transitions. The oscillator strengths of 
bound-bo und transition s consider ed here are from 
vALifllK upka et al.1 boool ). Kuruc^ JKurucz fc Belli Il995l ) 
and NIST on line databases. We take the vald log gf 
values as default for the whole transitions of an atom 
and choose the most acc urate one among vald, nist, 
iHirata fc Horaguchil ([19951), and Kurucz for the selected 
lines. For the bound-free transitions, we use the theoret- 
ical tables of photoionization cross-sections provided by 
the TopBase. We shift the TopBase energy threshold of 
photoionization to the corresponding nist values to keep 
consistency with the measured nist energy levels. In order 
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Figure 1. Photoionization cross-section as a function of the en- 
ergy of tlie incident photon for the 3s^ ^S ground state of Mgl. 
We use a slide window to smooth data after the resonance peak 
and undersample them. Dashed hne stands for data from Top- 
Base, red full line stands for smoothed and undersampled data 
used in our atomic model. 



to save computational time, we undersampled the pho- 
toionization cross-section tables taking care to keep large 
resonance peaks near the photoionization threshold but 
smoothing faint peaks far from the threshold using a slide 
window (see Fig. [T]). We assume that the photoionization 
cross-sections are the same for the fine structure of a level. 
Some photoionization data are missing for high levels, 
therefore we calculated photoionization threshold using 
hydro genic cross-section formulation of iKarzas fc Latten 
(|l96ll ). allowing us to compute Gaunt's factor as functions 
of the quantum numbers n and I. This approach leads 
to a better estimation than the sta n dard formulation of 
iMenzel fc Pekerid (|l935l ) and lMihalad (| 19781 ). 

CoUisional transitions. The inelastic collisions between 
the considered atoms with the other particles of the medium 
tend to set the atom in the Boltzmann Statistical Equilib- 
rium (SE). A careful treatment of these coUisional tran- 
sitions is therefore important for the SE that defines the 
level populations. Tables of quantum calculations of the col- 
lisional cross-sections with electrons are implemented when 
available. Otherwise, we use the mod ified Impact Param- 
eter Method (IPM) of ISeatonI l|l962al ') for allowed bound- 
bound collisions with electrons, which uses values of oscilla- 
tor strengths corresponding to radiative ones, and depends 
on the electronic d ensity and tem perature structures of the 
model atmosphere. ISeatonI (|l962al ) emphasized that the col- 
lisional rates are more reliable when (13.6/Ai5)/ is large, 
where AE is the energy of the transition in eV. Then, for 
intercombination lines with low / values, the coUisional rate 
can be under-estimated. For the forbidden transitions, when 
no quantum computations are available, it is generally as- 
sumed that the coUisional strength Q,ij does not depend on 
th e energy of the colliding electron s. We choose Q.ij = 1 (as 
in lAllenlll973l and lZhao et al.lll998l ) since quantum calcula- 
tions give this order of rnagnit ude for some forbidden optical 
lines (e.g. lBurgess et al.lll995l ). For forbidden transitions be- 
tween fine levels, we choose Q,ij = 10 in order to couple them 
by collisions. C oUisional ioniza tio n with electron s are com- 
puted following ISeatonI (|l962bh in lMihalaj l| 19781 ') where the 



coUisional ionization rates are taken to be proportional to 
the photoionization cross-sections. 

It is known today that the inelastic coUisional cross- 
sections with neutral hydrogen a tom s are largely over- 
estimated by the Drawin' s formula |Dra win 196 g) when ap- 
plied to dwarf stars (see | Caccin et al.M1993i . iBarklem et al.l 
|2003| . iBarklem et al.ll2011f ). This over-estimation tends to 
force LTE by exaggerating the coUisional rates with respect 
to the radiative ones. For giant and supergiant stars that 
we consider in this paper, we assume, a priori^ that the low 
gas density in the atmospheres of theses stars make the H- 
coUisions inefScient. Note that for the Mgi and Cai atoms 
efforts are being made to compute with quantum mechanics 
the exact cross-sections for in elastic collisions with neutral 
hydrogen (JGuitou et al.ll2010f ). These results will be imple- 
mented in a dedicated forthcoming paper. 

Damping parameters. The radiative damping is cal- 
culated using the Einstein coefficients of spontaneous de- 
excitation of each level present in the model atom: 



7rad = y^^ 



jk 



k<i 



^ik 



(1) 



where i stands for the lower level and j for the upper level. 
For the coUisional damping 7coii, we only take into account 
elastic collisions with neutral hydrogen, namely the Van 
der Waals' damping. Stark broadening is ignored in the 
present calculations, therefore: 7coii = 7vdw- We use quan- 
tum calculat ions of hydrogen e l astic c o Uisions with neutral 
species from'Anstce & O'Mara' ("iggsy jEaxklem fc O'Maral 
(1997); Barklom, O'Mara fc Ross (1991), hereafter named 
ABO theory. The great improvement of quantum mechani- 
cal calculation compared with the traditional Unsold recipe 
lUnsoldll 19551 ) is to lead to line widths that agree well with 
the ones observed without the need of an adjustable factor 
Fh (e.g. iMiha las 1978; Thcvcnin 1989 ). As used later in Ta- 
ble [3] and |4l we translate results from the ABO theory in 
terms of an enhancement factor Fh which is the ratio be- 
tween ABO broadening an classical Unsold broadening. We 
notice that for transitions including terms higher than I J5 4 
no broadening parameters are available from ABO theory. 
For these lines, we use an averaged value of this enhance- 
ment factor Fh ~ 2. We tested several values in the range of 
1 to 4 and found no significant effect on the resulting NLTE 
over LTE equivalent width ratios. 

2.2 Magnesium model atoms 

We developed two models for the neutral Mg atom: 

• one for the lines outside the Gaia/RVS wavelength 
range (model A) with the detailed fine structure of the en- 
ergy levels; 

• one dedicated to the multiplets in the Gaia/RVS wave- 
length range (model B) which need a special treatment to 
account for the blended multiplets rather than the separated 
components. 



2.2.1 Mgi model atom A 

Since the neutral Mg is an alkaline earth metal, its elec- 
tronic configuration leads to a singlet and a triplet systems. 
We consider all levels until n = IQ and I — 9 with fine 
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Figure 2. Mgl Grotrian diagram. The electronic configuration of the ground state is [Nc]3s^. Two multiplicity systems are present due 
to spin combination of the incomplete sub-level n = 3. All levels with fine structure are complete until n = 10 and I = A (but not visible at 
this scale); fine levels and multiplicity systems were merged for 5 ^ Z ^ 9. Only one doubly-excited state 3p^ is present in triplet system 
since in the singlet system the value of the corresponding energy level is greater than the ionization stage of Mgl. Horizontally dotted 
line stands for the ionization stage. Transitions represented by solid lines are those for which we deliver NLTE corrections; transitions 
represented by dashed lines are those in the Gaia/RVS wavelength range. 



structure for all levels with I ^ 4. For levels with I > 4, 
the multiplicity systems are merged (for the H, I, K, L and 
M terms) and no fine structure is considered. For complete- 
ness and because some high energy levels are missing in the 
database s, we used the f o rmula tion of polarization theory 
given by 'Chang fc Noved (|l983l ) to include them. We take 
into account 38 spectral terms in the singlet system and 38 
in the triplet system plus 15 merged terms (singlet -I- triplet 
systems without fine structure). In total, 149 energy levels 
+ ionization stage (at 7.646 eV) are considered for the SE 
computation (see the Grotrian diagram in Fig. [2]). 

1102 radiative bound-bound transitions have been 
taken into account: 127 allowed in the singlet system, 674 
in the triplet system and 301 intercombination transitions 
(semi-allowed transitions betwe en the two syste ms). We 
used hydrogenic approximation (jGreen et al.ll 19571 ) to com- 
pute several log gf missing in vald. For transitions be- 
tween levels with and without explicit fine structure, we did 
two assumptions to estimate individual oscillator strengths: 
firstly, as the transitions are produced between Rydberg 
states, the wavelength shifts between components are small 
enough to be neglected; and secondly the oscillator strength 
of each component of one multiplet is given in the LS 
coup hng. The TopBase photoionization cross-sections come 
from lButler et al.l (| 19931 '). 11026 electronic collisions between 
all levels and sub-levels are included using the IPM approx- 
imation when corresponding radiative transitions exist. The 
elastic coUisional broadening from the ABO theory is inter- 



polated for 87 transitions {s — p, p — d and d — /) in singlet 
and triplet systems, the Unsold recipe is used with Fh = 2 
for the remaining transitions. 



2.2.2 Mgi model atom B for the Gaia/RVS lines 

In order to reproduce Mgi lines in Gaia/RVS wavelength 
range, we construct a second Mgi model atom which en- 
ables MULTi to compute these lines as blended multiplets 
rather than separate components. These lines appear in the 
triplet system. We merge fine levels Ei under interest taking 
average energy level fi'ievei with respect to their statistical 
weights Qi: 



Ei, 



Y^=.9.Ef 



E»=i 9^ 



(2) 



as illustrated in Fig. |8] Averaged energy of each considered 
level and its fine structure are represented in Fig. [8] Then 
only averaged levels oi Ap ^P° , M ^D, Id ^D and 7/ ^F° 
are considered in this model atom B. The radiative transfer 
for a multiplet is then done under assuming that the same 
source function is used for the components of this multiplet 
(Mihalas 1978). Moreover we take the oscillator strength of 
the multiplet as a linear combination of the weighted gf of 
each component (in LS coupling). Actually, we can neglect 
the wavelength deviation between components and multiplet 
(relative deviation smaller than 1% for the Mgi lines): 
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Figure 3. Cal Grotrian diagram. The electronic configuration of the ground state is [Ar]4s^. All levels with fine structure are complete 
until n = 10 and 1 = 4. Twelve doubly excited terms are present: 4p^ ^SD ^P, 3d? ^PF (the terms for the singlet system ^SDG are 
missing but thought to be above the ionization threshold), SdAp^DF ^PDF° (the terms for the singlet system ^P° are missing) and 
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NLTE corrections; transitions represented by dashed lines are those in the Gaia/RVS wavelength range. 
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X 



(3) 



with respect to the line selection rules and with Xij (named 
Ritz wavelength) the average wavelength of the multiplet, 
gi = ^j gi, I for lower fine levels of term i, and u for upper 
fine levels of term j. We note that merging these fine levels 
implies shifts in wavelengths for the other lines. Then this 
atom will be exclusively used for the study of the lines in 
the Gaia/RVS wavelength range, except for the singlet line 
at 8473 A for which we use the Mg i model atom A. 



2.3 Calcium model atoms 

The neutral and first ionized Ca model atoms are built in 
the same way as the Mgi model atom A, but we consider 
them separately. We build sep arate model atoms of Ca i and 
Can, but for comparison with lMashonkina et al.l (l2007|) we 
merge the two separate atoms of Cai and Can to have a 
model atom of Ca through two ionization stages. 



Cai. This leads to consider 151 levels + ionization stage (at 
6.113 eV). A Grotrian diagram is shown in Fig. [3] Twelve 
doubly excited levels are presented for 4p^, 3d^, 3d4p and 
3d5s electronic configurations (see the legend of Fig. [3] for 
more detai ls). The TopBaso photoionization cross-sections 
come from ISaraph fc Storey (unpublished). All transitions 
present in vald are considered representing 2120 transitions: 
191 in the singlet system, 1161 in the triplet system and 768 
between these two systems. For elastic collisions with neu- 
tral hydrogen, the ABO theory provides broadening for 107 
lines. 

11628 electronic collisions are involved in the Cai model 
atom but only 11 collisional cross-sections with electrons for 
reson ance transitions are available from quantum computa- 
tions (jSamson fc Berringtonl l2001). These authors used the 
Variable Phase Method (VPM) for Cai lines shown in Ta- 
ble O For inelastic collisio ns with neutral hydrogen, we can 
refer to lMashonkina et al.l (120071 ) who show detailed tests of 
the importance of these inelastic collisions. They conclude 
that for many transition lines these collisions are small for 
Cai and negligible for Can to reproduce the line profiles. 



2.3.1 Cai model atom 

Neutral Ca is also an alkaline earth metal with two multiplic- 
ity systems: singlet and triplet systems. However, more dou- 
bly excited energy levels exist which produce a more complex 
model atom than the Mg i one. All energy levels until n = 10 
and I = 4 are considered with fine structure and we take into 
account all doubly excited levels under ionization stage of 



2.3.2 Can m,odel atom 

The Can atom is like an alkali element, so it presents 
only one system of levels with a doubly multiplicity sys- 
tem (see FigUJ. We consider levels until n = 10 and / = 5 
for Can. Fine structure is used, leading to 74 levels for 
Can. 422 radiative transitions are considered from the vald 
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Figure 4. Call Grotrian diagram. First ionized calcium has a 
doubly multiplicity system. In this case, we have represented all 
the radiative bound-bound transitions used in the model atom 
with full lines. Transition represented by red lines are those for 
which we will deliver NLTE corrections. Horizontally dotted line 
represents the ionization stage. 



Ap^P° of Ca n. This model atom is use d to compare results 
with those of iMashonkina et al.l (|2007l ') in Sect. 14.4.11 



3 TEST ON TWO BENCHMARK STARS: THE 
SUN AND ARCTURUS 

The new model atoms of these a-elements are first tested 
on two well known benchmark stars: a main sequence star 
(the Sun) and a giant (Arcturus). We compute flux profiles 
for several lines in the Gaia/RVS wavelen gth range usin g 
the NLTE radiative transfer code multi ( CarlssonI Il98q l. 
This code directly provides correction ratios of NLTE {W) 
versus LTE {W*) equivalent widths (EW), hereafter named 
NLTE/LTE EW ratios or W/W*. 

We define the abundance of an eleme nt (El) relative 
t o hy drogen by Aei using notatiorjj from iGrevesse et al.l 
( 20071 ). Stellar abundances are given with respect to the so- 
lar ones on a logarithmic scale [El/H] as usually donqj. The 
abundances of a-elements are given with respect to iron by 
[a/Fe]. 

We also define NLTE abundance correction as: 



A [El/H] = ylii^™ 



(4) 



where ^4^1 is the LTE abundance of the element El, and 
A^Y^'^ the NLTE abundance. 

For very weak lines (e.g. Mgi 8473 A) or for very strong 
lines (e.g. CaT), one can express W/W* as a function of the 
abundance ratio by using the curve of growth theory: 



W 
W* 



A'" 



4 LTE 
^El 



(5) 



where /3 is a coefficient which depends on the strength of 
the line. If log(VK/A) ^ —4.8, the line is considered weak 
and \og{W/X) has a linear dependence on the abundance 
(/3 = 1); if log(VF/A) ^ —4.2 the line is considered strong and 
log(M^/A) has a square root dependence on the abundance 
(/3 — 1/2). We deduce for these two extreme cases the NLTE 
abundance correction as: 



database. Photoionization cros s-sections from the TopBase 
ijSaraph fc StorevI lunpublishedh do not present any peaks. 
In the ABO theory, the line broadening due to elastic col- 
lisions with_Ji^djogen_Js_onW_known for the H, K and CaT 
lines l|Barklem fc O'Maral Il998l l. 2775 inelastic collisions 
with electrons are involved in the Ca n model atom . Quan- 
tum computations are available from lBurgess et al.l ([1995|), 
using the Coulomb Distorted Wave approximation (CDW), 
for 21 transitions between mean levels: 10 allowed and 11 for- 
bidden as reported in Table O For components between fine 
levels, we apply the same ratio from the oscillator strengths 
gf to the coUisional strength Q,. 



2.3.3 Ca i/ii model atom 

In order to see if NLTE effects act on the ionization equilib- 
rium, we build a Ca i/ii model atom by merging the Ca i and 
the Can ones. Several energy levels of Cai are photoionized 
on excited levels of Can. It is the case for the doubly ex- 
cited levels 3d4p^DF ^PDF° which photoionize on to Sd^D 
of Can and for the 4p^ ^SD "^P which photoionize on to 



A [El/H] = --log(M//W* 



(6) 



We apply this relation to deduce NLTE abundance correc- 
tions A [EL/H] from W/W* when we compare weak and 
strong lines with the available literature as done in Table [U 
[6] and [9l Note that between these two extreme cases, /? is 
lower than 1/2. This can imply strong variations of A [El/H] 
corresponding to the plateau of the curve of growth (e.g. 
iMihalad [1973 . pp 319, 320 and 325). We emphasize that 
Eq. [6] is an approximation that gives an estimate of the 
NLTE abundance correction. The right way to determine 
the NLTE abundance correction should be to vary the NLTE 
abundance of El until the NLTE EW fits the observed EW. 
In the following sections, we test our model atoms on 
the Su n for which we know in advance that NLTE effects are 
small (|Asplundll2005l ) but that can be used to estimate the 
missing enhancement factor Fh for some Gaia lines (outside 
of the ABO theory range). We also test our model atoms on 



^ ^El = logiQ (Nei/Nh) -f 12, N the number density of a given 
element. 

6 [El/H] = A*i - A%. 



A grid of NLTE corrections for Mg and Ca for evolved stars 7 




0.2 




— 


Solar FTS 

MULTI LTE (MARCS) 
MULTI NLTE (MARCS) 


87 


I)A 


8735 8736 
A [A] 



0.4 1 L 

8708 




Solar FTS 

MULTI LTE (MARCS) 

MULTI KLTE (MARCS) 

J \ , , , \ , , , \ , , , \ 

8712 8714 8716 8718 





7 [A] 




1.2 


- 


- 


1.0 


^— ^.^^^ 


__ — —7 


0.8 


: ^%\ /^ 


y - 


0.6 


: V\ / 


- 


0.4 


: y 


^ 


0.2 


- Solar ETS " 

MUl Tl 1 TF fMARCSl 


^ 




MUl Tl MITE fMARCSl 


- 


00 


1 , , , 1 , , 


, 



0.0^ 
8734 



0.4 I L 

8703 




Arcturus 

MULTI LTE (MARCS) 

MULTI NLTE (MARCS) 



... Solar ETS 

_ MULTI LTE (MARCS) 

_ MULTI NLTE (MARCS) 
J \ , , , \ , , , \ , 

8712 8714 8716 




Solar FTS 

MULTI LTE (MARCS) 

MULTI NLTE (MARCS)' 



8498 
A [A] 



0.0 I 

8494 



8498 

>' [A] 



Figure 5. Solar line fits in the Gaia/RVS wavelength range, using 
a solar plane— parallel marcs model. The blue lines stand for LTE 
synthesis and the red lines for NLTE synthesis. Dashed lines stand 
for observations. Top panel: Mgl 8736 A multiplet. Middle panel: 
Mgl triplet at AA 8710, 8712 and 8717. Bottom panel: Call 8498 
A line. 



Arcturus, a well known red giant, using high resolution and 
high S/N spectrum. 



Figure 6. Arcturus line fits in the Gaia/RVS wavelength range, 
using a spherical interpolated MARCS model atmosphere. The 
blue lines stand for LTE synthesis and the red lines for NLTE 
synthesis. Dotted lines stand for observations. Top panel: Mgl 
8736 A multiplet. Middle panel: Mgl triplet at AA 8710, 8712 
and 8717. Note that NLTE effects are more pronounced than the 
Sun (Fig. [Sjl. Bottom panel: Can 8498 A line. Departure from 
LTE appears in the core. 



3.1 NLTE synthetic line profile fits of solar lines 

We compare our NLTE disk integrated line profiles with 
the observed solar spectrum obtained by the Fourier Trans- 



8 T. Merle, F. Thevenin, B. Pichon, L. Bigot 



Table 1. Solar NLTE abundance corrections for Mgl lines of astrophysical interest with (A[Mg/H]//, Sh = 1) and without (A[Mg/H], 
Sff = 0) inelastic collisions with neutral hydrogen according to a scaled Drawin's formula. Wc applied Eq.lJGjl to compute A[Mg/H] and 
A[Mg/U]H. 





W [mA] 


Wh [mA] 


W/W 


[mA] 


w/w;^ 


This work 




Zhao et al. 


(1998) 


A [A] 


A[Mg/Hj 


AlMg/H|« 


A[Mg/Hj A 


[Mg/Hjn 


4571.095 


83 


85 




0.935 


0.954 


0.06 




0.02 


0.07 


0.03 


5167.320 


675 


706 




0.954 


0.998 


0.04 




0.00 


0.04 


0.01 


5172.683 


1139 


1195 




0.951 


0.998 


0.04 




0.00 


0.04 


0.01 


5183.603 


1430 


1502 




0.949 


0.998 


0.05 




0.00 


0.04 


0.01 


5528.403 


333 


342 




0.973 


1.001 


0.02 




0.00 


0.06 


0.00 


8806.751 


465 


499 




0.934 


1.001 


0.06 




0.00 


0.09 


-0.03 


8923.563 


43 


47 




0.929 


1.001 


0.03 




0.00 


0.11 


0.02 


11828.18 


755 


771 




0.981 


1.001 


0.02 




0.00 


0.05 


0.02 



form Spectrograph, here-after FTS (Brault & Neckel 19870 
This is particularly well adapted for the present work be- 
cause of its high signal-to-noise ratio {S/N > 500) and its 
large resolving power {R ~ 400000). We use a ID, static, 
and plane-parallel solar model atmosphere from marcs grid 
(Gustafsson et al. 2008) with T^ff = 5777 K a nd log g = 4.44 
The so lar abundances of Ca and Mg are from lGrevesse et al.l 
l|2007h corresponding to A®^ = 7.53 and A®^ = 6.31. The 
flux profiles are convolved with a gaussian profile to account 
for the effects of macroturbulence velocity (~ 2 km s~^). 
As the projected rotational velocity is weak for the Sun, we 
neglect it. 

The NLTE effects on lines in the visible wavelength 
range for both Mgi and Cai are small in the late type 
solar-like stars (e.g . comparing the NLTE abundances from 



Idiart fc Theveninl i2000 with the LTE ones from Theveni n 
19981 ). However, NLTE effe cts become non-n egligible for IR 



lines (see the conclusion of IZhao et al.lll998l ). As an exam- 
ple, we compare NLTE and LTE profiles for Mgi lines in 
the RVS wavelength range. In the top panel of Fig. [SI we 
show the Mgi 8736 A multiplet computed with the model 
B. This line is fitted in NLTE with a value of Fh = 3 and 
^Mg™ = 7.63. We find a ratio of W/W = 1.10 and a NLTE 
EW of 204 mA. The LTE fit gives a value of Al:J^ = 7.70. 
Therefore, we have A[Mg/H] = -0.07 dex of NLTE abun- 
dance correction for this line on the Sun. The middle panel 
of Fig. [S] presents the Mgi triplet at AA 8710, 8712 and 8717 
computed assuming the same source function since we have 
merged lower and upper fine levels into one lower and one 
upper in model B. Our best fit is obtained with Fh = 2.5, 
^Mg^"^ = 7.63 and we found that NLTE effects are negligible 
for this multiplet since W/W* ^ 1.02. 

We compare our NLTE results for the Sun with 
I Zhao et al.l l|l998h wh o used Kurucz's model atmosphere 
with the DETAIL code (jButler fc Giddine:slll985l ) to perform 
NLTE computations. They studied the impact of the in- 
elastic collision with hydrogen on the NLTE results, using 
Drawin's formula and empirical scaling laws. We compute 
A[Mg/H] using Eq. ((5]) for several lines in common with 
their study. For the comparison, we add collisions with neu- 
tral hydrogen, using Drawin's formula with a scaling factor 
Sh ~ 1. Results are presented in Table [T] for solar lines. 
For the Mgi b t riplet we obtain a very good agreement 
with IZhao et all (| 19981 '). For other transitions, our NLTE 
corrections are lower than theirs. The effect of inelastic col- 



lisions with neutral hydrogen, through Drawin's formula, 
eliminates completely the NLTE effects except for the in- 
tercombination resonance li ne at 4571 A for w hich we have 
A[Mg/U]H = 0.02 (0.03 in lZhaoet ai]|l998l ). In the Sun, 
NLTE corrections for these lines are all less than 0.06 dex. 
If we do not take into account Drawin's formula, all the 
NLTE corrections are positive. We have already mentioned 
in Sect. 12.11 that Drawin's formula is not adapted for the 
inelastic collisions with neutral hydrogen. 

We also consider NLTE effects on the well studied CaT. 
We note that the central depression of these lines are deeper 
in NLTE compared to LTE ones as seen in the bottom panel 
of Fig. [5] The discrepancies with the observations are due 
to the fact that the cores of these lines form in the chro- 
mosphere which is not included in the marcs model atmo- 
spheres. Since these strong lines are dominated by the wings 
that are formed in LTE, the equivalent widths are quite 
insensitive to NLTE in solar conditions {W/W* ~ 1.00). 
Note that the small asymmetry of th e line c o re ma y be 
due to an isotopi c shift as shown by iKuruca (2005) and 
iMashonkina et al.l ((20071) but a contribution of 3D effect is 
also not to be excluded. 



3.2 NLTE synthetic line profile fits of Arcturus 

We compute NLTE line profiles using an interpolated marcs 
model atmosphere (spherical geometry, Tcef = 4300 K, 
logg = 1.5, [Fe/H] = -0.5, [o/Fe] = +0.2, ^ = 2 km s'S 
M = \ Mq) obtained with the interpolation code of 
T. Masserorjj. The spectrum of Arcturus comes from the 
Narval instrument at the Bernard Lyot Telescopqj {R ~ 
65000) with hight signal-to-noise ratio {S/N > 500). 

The computation of Arcturus line profiles reveals that 
NLTE effects are larger than in the Sun in particular for the 
Gaia/RVS lines as shown in Fig. (6] When we fit the Mg i 
8736 A line with A^l^^ = 6.93, we find an equivalent width 
of 128 mA and a ratio of W/W* = 1.30. The LTE fit of the 



same line gives a value of A 



Mg 



7.23. Therefore, we have 



A[Mg/H] = -0.30 dex of NLTE abundance correction for 
this fine in Arcturus. For the Mgi triplet at AA 8710, 8712 
and 8717, the NLTE effects are smaller than for 8736 A line 
but larger than in the Sun. We find W/W* = 1.10. 

Concerning the Can IR 8498 A line, departure from 



* available at: 
'^ available at: ftp://ftp.hs. uni-hamburg.de/pub/outgoing/FTS-Atlas/| ^ available at: 



http://marcs.astro.uu.se/software.phpl 



http://magics.bagn.obs-mip.fr| 
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LTE always affects the line core. We are not able to repro- 
duce satisfactorily the observed line core, neither in LTE or 
NLTE as shown in bottom panel of Fig. |6] Contrary to the 
solar case, the central depression of these lines are deeper in 
LTE compared to NLTE ones. The best value of Ca abun- 
dance fitting the wings of this hue in NLTE is ^ca^"^ = 5.81. 
The equivalent width ratio is W/W* = 0.95. As in the Sun, 
the discrepancy in the line core is probably due to the pres- 
ence of a chromosphere in such giant stars and in particular 
in Arcturus for w hich many studies have been devoted (e.g. 
lAvres et al]| 19821 ). 



4 NLTE EFFECTS VERSUS ATMOSPHERIC 
PARAMETERS 

We compute an NLTE correction table for a wide range 
of stellar parameters (Tcff, logg, [Fe/H]) of cool giant and 
supergiant stars with metallicities from -1-0.5 down to —4 
dex. The selected lines are those of astrophysical importance 
in the optical and near IR wavelength ranges, especially for 
those in the domain of the Gaia/RVS. Then, we describe 
the general trends of NLTE effects for the selected Mg i, Ca i 
and Ca ii lines of astrophysical interest (Table [3] and |4)| as 
functions of stellar parameters. The most important results 
deduced from our computed W/W* are an enhancement of 
NLTE effects with decreasing metallicity, a non monotonic 
variation with the effective temperature, and a significative 
dependence on the surface gravity for several lines. We also 
compare our results with previous works. 



4.1 The W/W* grid for giant and supergiant stars 

We compute NLTE/LTE EW ratios for T^ff , log g and [Fe/H] 
covering most of the late-type giant and supergiant stars. 
The effective temperature Tos range is 3500(200)3900 K 
and 4000(250)5250 K, the surface gravity \ogg range 
is 0.5(0.5)2.0 dex, and the metallicity [Fe/H] range is 
-4.00(0.50) - 1.00(0.25) + 0.50 dex. We use spherical 
ID, and static model a tmospheres from the marcs sitqij 
IjCustafsson et al.l 20081) with c hemical standard composi- 
tion from ICrevesse et al.l (|2007l ). In the standard composi- 
tion of MARCS models, [a/Fe] depends on [Fe/H] as: 



[Fe/H] > 


^ [a/Fe] = 







[Fe/H] < 


^ [a/Fe] = 


-0.4 [Fe/H] 


(7) 


[Fe/H] < 


-1 ^ [a/Fe] = 


-hO.4 





-1 < 



This means that the Ca abundance used for NLTE compu- 
tations follows the [a/Fe] value of the model atmosphere. 
The MARCS models with M = 1 Mq and constant microtur- 
bulent velocity fields ^ = 2 km s~^ are adopted. We inter- 
polated MARCS model atmospheres with the Masseron code 
for [Fe/H] = —3.50 and for few models missing in the marcs 
database. 

For the NLTE computations, we need to solve coher- 
ently the coupled equations of the SE and the radiative 
transfer which use the local approximate lambda operator 
implemented in the ID, plane -parallel, and t ime indepen- 
dent code MULTi, version 2.2 IjCarlssonl Il986l ). Continuous 



Table 2. Part of the results for the Mgl 8736 A and Call 8498 
A lines. The full version of the results is available in the online 
version of the journal. NLTE/LTE EW ratios W/W* are given 
for Toff = 5000 K, logg = 1 and 2, and metallicities if VK > 1 mA. 
We do not take into account collisions with neutral hydrogen 
(Sjj = 0). The NLTE abundance correction A[E1/H] is given 
when the Eq.lJGJI is applicable. 



T,ff[K] logs [Fe/H] [a/Fe] 



Mgl 8736.012 A 

W[mA] W/W* A[Mg/H] 



5000 



+0.50 


0.0 


225 


1.179 




+0.25 


0.0 


205 


1.226 




+0.00 


0.0 


177 


1.268 




-0.25 


0.1 


158 


1.299 




-0.50 


0.2 


135 


1.302 


-0.11 


-0.75 


0.3 


111 


1.272 


-0.10 


-1.00 


0.4 


89 


1.210 


-0.08 


-1.50 


0.4 


33 


0.927 


+0.03 


-2.00 


0.4 


9 


0.667 


+0.18 


+0.50 


0.0 


211 


1.170 




+0.25 


0.0 


192 


1.196 




+0.00 


0.0 


167 


1.214 




-0.25 


0.1 


151 


1.230 




-0.50 


0.2 


131 


1.236 


-0.09 


-0.75 


0.3 


110 


1.218 


-0.09 


-1.00 


0.4 


88 


1.170 


-0.07 


-1.50 


0.4 


34 


0.919 


+0.04 


-2.00 


0.4 


10 


0.683 


+0.17 



T,ff[K] logg [Fe/H] [a/Fe] 



Can 8498.018 A 
iy[mA] W/W* A[Mg/H] 



5000 



+0.50 


0.0 


3857 


0.969 


+0.03 


+0.00 


0.0 


3086 


0.973 


+0.02 


-0.50 


0.2 


2550 


0.983 


+0.01 


-1.00 


0.4 


2036 


0.993 


+0.01 


-2.00 


0.4 


863 


1.037 


-0.03 


-3.00 


0.4 


396 


1.154 




-4.00 


0.4 


264 


1.352 




+0.50 


0.0 


2255 


0.965 


+0.03 


+0.00 


0.0 


1909 


0.972 


+0.02 


-0.50 


0.2 


1715 


0.984 


+0.01 


-1.00 


0.4 


1446 


0.998 


+0.00 


-2.00 


0.4 


737 


1.069 


-0.06 


-3.00 


0.4 


370 


1.233 




-4.00 


0.4 


241 


1.435 





available at: |http://niarcs. astro. uu.se| 



opacities come from Uppsala package (JGustafssonl Il973l l . 
We take into account 45000 lines for line opacity compu- 
tations. Opacities from molecular lines are not included. We 
do not take into account the collisions with neutral hydrogen 
{Sh = 0). Starting population numbers are taken at LTE 
and are used to be compared with the NLTE results. With 
very large model atoms, the convergence of the Newton- 
Raphson scheme can be very slow and the number of iter- 
ations for some extremely metal-poor models can be very 
high. For the Mg i and Ca ii model atoms, the minimal rela- 
tive precision for the population number at any layer of the 
atmosphere is of 10"* whereas it decreases to 10"'^ for Cai 
for which the model atom is more complex. 

We emphasise that we did ID NLTE plane-parallel 
radiative transfer in ID theoretical spherical model atmo- 
spheres which may appear not consistent. However, as shown 
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by iHeiter fc Eriksson! l|2006l '). it is better to use spheri- 
cal model atmosphere with plane-parallel radiative transfer 
rather than plane-parallel model atmosphere with plane- 
parallel radiative transfer for giant and supergiant stars. 

Moreover the W/W* do not take into account the pos- 
sible contribution of the stellar chromospheres for lines such 
as the Mgi b triplet, the Can H, K and the CaT. AU the 
MARCS models are theoretical, ID and LTE models with- 
out any model of chromosphere, i.e. without any rising of 
the temperature profile in the upper layers. This absence 
of model chromosphere can explain why we are not able to 
reproduce the line core of CaT in Arcturus (Fig. [S]). For an 
extended discussion about the treatment of chromospheric 
contribut ion on a peculiar r ed giant star (/3 Get), using this 
work, see lBerio et al.1 (|201ir ). 

The NLTE results for the lines in Tables [3] and g] are 
available in electronic forms. Illustrations of W/W* for a 
selection of these lines as function of [Fe/H] , for several ToEF 
and parametrized in loggr are shown in Appendix A and B. 
Global NLTE results are summarized in Tables[5l[7]and[8]for 
Mgi, Gai and Gaii lines, respectively. We show in Table [21 
an example of the NLTE results expressed in terms of W/W* 
for the Mgi 8736 A line and the Gaii IR 8498 A line as a 
function of stellar parameters. We also give the NLTE EW in 
order to know if the Eq. (|6]) is applicable to directly deduce 
the NLTE abundance correction A[E1/H]. 

4.2 Line selection 

The lines are mostly selected as unblended and strong 
enough to be measurable in metal-poor stars. Spectroscopic 
and micro-physics details on Mgi, Gai and Gaii selected 
lines are presented in Table [31 

Wavelengths are given in air using the lAU standard 
dispersion law for conversion between vacuum and air (see 
lMortonlll99ll ). For Mgi, we mainly selected lines of the sin- 
glet system except for the Mgi b triplet lines at AA 5167, 
5172, 5183, at 7657 and 8997 A lines, and the lines belonging 
to the Gaia/RVS domain. The important intercombination 
resonance line at 4571 A is also selected. This line has been 
recently studied ([Langangen fc Garlssonll2009l ') and has been 
proposed to be used as temperature diagnostic in stellar at- 
mospheres. For Gai, we mainly selected lines coming from 
singly excited levels except for the 5512 A line and for the 
weak triplet in RVS because for lines coming from doubly 
excited levels we have less atomic data. Most of these opti- 
cal lines are blended in their wings at solar metallicity but 
are less blended with decreasing metallicity. For the optical 
red triplet at AA 6102, 612 2, 6162, we p r efer t o use the en- 
hancement factors Fh from lGavrel et al.l (|l99y) who provide 
resolved Fh for each transition of the triplet {Fh = 2.92, 
2.22 and 2.44 respectively) rather than ABO theory which 
gives one unique value of Fh ~ 2.55. Lines selected for Gaii 
are the strong H & K lines, the IR triplet and few other 
unblended IR lines. 

Lines in the Gaia/RVS are dominated by the GaT. 
There are also 5 Mgi lines: a very weak line at 8473 A, 
a weak triplet at AA 8710, 8712, 8717, and a stronger line 
at 8736 A. Details on the components of these lines are 
shown in Table [31 A partial Grotrian diagram in Fig. [51 
also details the fine structure of the levels involved in these 
multiplets. Note that only the mean level of the 7/ "^F" 
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Figure 7. Mg I departure coefficient as a function of optical depth 
in the continuum at 5000 A for a metal-poor giant stellar model. 
We separate levels as function of their energy Ei with respect to 
the ground state as shown in the legend. 



is available on the nist database. The weakest 8710 and 
8712 A components of the triplet are also blended on their 
red wings by iron lines as shown in Fig. [5l and [6l This mul- 
tiplet at AA 8710, 8712, 8717 is represented, using the Mgi 
B model, with a Ritz wavelength at 8715 A for which we 
deliver NLTE/LTE EW ratio by assuming the same NLTE 
behaviour for the three components. For Gai, 5 transitions 
exist in this domain, between doubled excited levels 3d4p 
and 3d5s, as shown in Fig. [31 These lines are very weak 
at solar metallicity and only AA 8525, 8583 and 8633 are 
visible in the solar spectrum. Moreover, in the spectrum 
of Arcturus, these three lines are a little bit strengthened 
but blended by CN lines. Thus, we decide to provide NLTE 
corrections only for these three lines. We notice that for 
the AA 8498, 8542, 8662 IR triplet (CaT), we used theoreti- 
cal log q.f values ( — 1.3 56, —0.405, —0.668 respectively) from 
iMelendez et al.l ((2007]) which are in a very good agreement 
with values derived from 3D hydrodyn amical line fits on 
the Sun (-1.309, -0.410, -0.683) from JBigot fc Thevenin! 
(,200& ). 



4.3 The Mgi lines 

In order to see how the Mgi is affected by the NLTE, the 
departure coefficients hi = nt/n* of the Mgi model atom are 
plotted in Fig. [3 where rii and n* stands for NLTE and LTE 
level populations respectively, as a function of the standard 
optical depth rsooo for a metal-poor giant model with atmo- 
spheric parameters of Teff = 4500 K, logg = 2, [Fe/H] = — 3 
and [a/Fe] = -1-0.4. The ground state is very depopulated 
regarding the LTE case. This is due to an increase of the 
UV radiative field produced by a decrease of the metallic- 
ity. Moreover, we show that there is an overpopulation of a 
large part of higher levels explained by a collisional domi- 
nated recombination. The population in the ionization stage 
is in LTE. 

The Mgi AA 4167, 4702, 5528 and 8806 lines are tran- 
sitions between the 3p ^P° level and ^D terms. Therefore, 
their W/W* have similar behaviours except for the 8806 A 
line which has the most pronounced deviation from LTE (see 



A grid of NLTE corrections for Mg and Ca for evolved stars 11 

Table 3. Mgl, Cal and Call line selection for NLTE corrections. The wavelengths A are given in the air. The oscillator strengths gf and 
the radiative damping parameters 7rad are given in a decimal logarithmic scale. The enhancement factor Fu is defined as the ratio of 
the coUisional damping parameter computed from ABO theory divided by the classical Unsold's formula (a default value of 2 is applied 
when lines are outside the range of ABO theory). The radiative damping parameters are calculated using the log gf values taken from 
VALD, except when specified. Mgl line selection in the Gaia/RVS wavelength range is detailed in Table |4] The precision of log gf are 
given with letters, following the NIST notation: A (< 3%), B+ (< 7%), B (< 10%), C+ (< 18%), C (< 25%), D (< 50%), and E ( > 50%). 
Inelastic collisions with electrons are specified in the last column (Coll.) where IPM stands fo r Impact Parameter Method from lSeatoiJ 
(ll962a|), CDW stan ds for Coulomb Distor t ed W ave approximation from lBurgess et al.l l|l995l ) for Can, and VPM stands for Variable 
Phase Method from lSamson &: BerringtonI 11200 Ih for Cal resonance lines. 



A [A] 



Lower level 
Configuration 



g^ E, [eV] 



Upper level 
Configuration 



9j Ej [eV] log gf Prec. log 7,,ad 



Fh 



Coll. 



Mgl 



4167.271 


3p 


\ po 


3 


4.345802 


Id 


Id 


5 


7.320153 


-0.745" 


C+ 


8.71 


2.00 


IPM 


4571.095 


4s^ 


is 


1 


0.000000 


3p 


3 po 


3 


2.711592 


-5.688'' 


C+ 


2.34 


1.34'» 


IPM 


4702.990 


3p 


1 po 


3 


4.345802 


bd 


Id 


5 


6.981349 


-0.440'" 


B+ 


8.72 


2,770, 


IPM 


4730.028 


3p 


1 po 


3 


4.345802 


6s 


is 


1 


6.966284 


-2.400'' 


B+ 


8.71 


2.00 


IPM 


5167.320 


3P 


3po 


1 


2.709105 


4s 


^s 


3 


5.107827 


-0.931* 


B 


7.98 


2.37'* 


IPM 


5172.683 


3p 


3po 


3 


2.711592 


4s 


3s 


3 


5.107827 


-0.450* 


B 


7.98 


2.37'* 


IPM 


5183.603 


3p 


Zpo 


5 


2.716640 


4s 


3s 


3 


5.107827 


-0.239* 


B 


7.98 


2.37" 


IPM 


5528.403 


3p 


1 po 


3 


4.345802 


Ad 


Id 


5 


6.587855 


-0.490'' 


B+ 


7.38 


2.11" 


IPM 


5711.086 


3p 


1 po 


3 


4.345802 


5s 


is 


1 


6.516139 


-1.720'' 


B+ 


8.71 


2.00 


IPM 


7657.599 


4s 


3s 


3 


5.107827 


5p 


3 po 


5 


6.726480 


-1.268" 


C 


8.00 


2.00 


IPM 


8473.688 


4p 


3po 


5 


5.932787 


9s 


3s 


3 


7.395551 


-2.020'= 




7.20 


2.00 


IPM 


8806.751 


3p 


1 po 


3 


4.345803 


3d 


ID 


5 


5.753246 


-0.134° 


A 


8.72 


1.39" 


IPM 


8923.563 


4s 


15 


1 


5.932370 


5p 


1 po 


3 


7.354590 


-1.678° 


B 


7.81 


2.00 


IPM 


8997.147 


4p 


Zpo 


5 


5.932787 


8s 


3S 


3 


7.310447 


-1.770 




7.24 


2.00 


IPM 


10312.52 


4p 


1 po 


3 


5.945915 


7d 


Id 


5 


7.364755 


-1.730° 


C+ 


8.19 


2.00 


IPM 


11828.18 


3p 


1 po 


3 


4.345802 


4s 


is 


1 


5.753246 


-0.333'" 


A 


8.73 


2.91" 


IPM 



Cai 



4226.727 


4s2 


is 


1 


0.000000 


4p 


1 po 


3 


2.932512 


0.244 


B+ 


8.34 


2.09" 


VPM 


4425.435 


4p 


3po 


1 


1.879340 


4d 


3d 


3 


4.680180 


-0.358 


C 


7.94 


1.78" 


IPM 


4578.549 


3d 


3d 


3 


2.521263 


4.f 


3^0 


5 


5.228440 


-0.697 


B+ 


7.36 


2.00 


IPM 


5512.980 


4p 


1 po 


3 


2.932512 


4p2 


is 


1 


5.180837 


-0.447'' 


B+ 


8.48 


2.00 


IPM 


5867.559 


4p 


1 po 


3 


2.932512 


6s 


is 


1 


5.044971 


-1.570 


B 


8.35 


2.00 


IPM 


6102.719 


4p 


3po 


1 


1.879340 


5s 


3s 


3 


3.910399 


-0.850* 


B 


7.87 


2.92*3 


IPM 


6122.213 


4p 


3po 


3 


1.885807 


5s 


3s 


3 


3.910399 


-0.380* 


B 


7.87 


2.22*3 


IPM 


6161.294 


3d 


3d 


5 


2.522986 


5p 


3 po 


5 


4.534736 


-1.266 


B+ 


7.30 


2.11" 


IPM 


6162.170 


4p 


3po 


5 


1.898935 


5s 


3s 


3 


3.910399 


-0.170* 


B 


7.87 


2.44*3 


IPM 


6166.438 


3d 


3d 


3 


2.521263 


5p 


3po 


1 


4.531335 


-1.142 


B+ 


7.28 


2.11" 


IPM 


6169.037 


3d 


3d 


5 


2.522986 


5p 


3 po 


3 


4.532211 


-0.797 


B+ 


7.30 


2.12" 


IPM 


6169.562 


3d 


3d 


7 


2.525682 


5p 


3 po 


5 


4.534736 


-0.478 


B+ 


7.30 


2.11" 


IPM 


6455.592 


3d 


3d 


5 


2.522986 


3d4p 


lp)o 


5 


4.443025 


-1.290'' 


B+ 


7.67 


2.00 


IPM 


6572.777 


4s2 


is 


1 


0.000000 


4p 


3 po 


3 


1.885807 


-4.296'' 


D+ 


3.47 


1.57" 


VPM 


7326.142 


4p 


1 po 


3 


2.932512 


4d 


Id 


5 


4.624398 


-0.208 


B 


8.37 


1.72" 


IPM 


8525.707 


3d4p 


3po 


5 


4.430011 


3d5s 


3d 


3 


5.883850 


-1.202'= 




7.70 


2.00 


IPM 


8583.318 


3d4p 


3po 


7 


4.440954 


3d5s 


3d 


5 


5.885035 


-0.986'= 




7.79 


2.00 


IPM 


8633.929 


3d4p 


3po 


9 


4.450647 


3d5s 


3d 


7 


5.886263 


-0.810'= 




7.77 


2.00 


IPM 



Can 



3933.663 
3968.468 
8248.791 
8498.018 
8542.086 
8662.135 
8912.059 
9854.752 
11949.74 



4s 
4s 
5p 
3d 
3d 
3d 
4d 
5p 
5s 



2p 
2p 
2po 

2d 

2d 

2d 

2d 
2po 

2S 



0.000000 
0.000000 
7.514840 
1.692408 
1.699932 
1.692408 
7.047168 
7.505137 
6.467874 



4p 
4p 
5d 
4p 
4p 
4p 

4/ 
6s 
5p 



2po 
2po 

2d 
2po 

2po 

2po 

2po 

2S 
2po 



3.150984 
3.123349 
9.017485 
3.150984 
3.150984 
3.123349 
8.437981 
8.762907 
7.505137 



0.104'' 
-0.210'' 

0.556 
-1.356'^ 
-0.405'^ 
-0.668'^ 

0.575'= 
-0.260° 
-0.010 



A4 
A^ 

c 



D 
D 



8.21 


1.42" 


CDW 


8.15 


1.43" 


CDW 


8.29 


2.00 


IPM 


8.21 


2.16" 


CDW 


8.21 


2.16" 


CDW 


8.15 


2.16" 


CDW 


8.85 


2.00 


CDW 


8.26 


2.00 


IPM 


8.43 


2.00 


CDW 



''iHirata fc Horaguchil lll995l '). 

set al.N200g|) for Cal 

ABO theory '^ ICavrel etld] l ll996l'l . 



° NIST, 

lAldenii 



■■ iKurucz fc Bell lll995l '). "^ iMelendez et al.l | |2007|) . * lAldenius et"aLl ll2007t) for Mg I and 
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Table 4. Details of the blended multiplets of Mgl in the Gaia/RVS wavelengths. The log gf are from VALD except when specified. The 
wavelengths in bold and between (.) stand respectively for observed lines and Ritz wavelengths (theoretical ones). The notations are the 
same as in Tabic [3] 





Lower level 






Upper level 
















A [A] 


Configuration 


9i 


E, [eV] 


Configuration 


9] 


£, [eV] 


logs/ 


Prec. 


log7r 


Fh 


Coll. 


(8715.255} 


4p 


3 po 


9 


5.932370 


7d 


3d 


15 


7.354590 








2.50 


IPM 


8710.169 


4p 


3 DO 


1 


5.931542 


7d 


3d 


3 


7.354592 


-1.565" 


D 


7.38 


2.50 


IPM 


8712.678 


4p 


3po 


3 


5.931951 


7d 


3d 


8 


7.354591 


-1.088 


D 


7.39 


2.50 


IPM 


8712.670 






3 


5.931951 






3 


7.354592 


-1.690" 


D 


7.39 






8712.684 






3 


5.931951 






5 


7.354590 


-1.213" 


D+ 


7.32 






8717.812 


4p 


3 po 


5 


5.932787 


7d 


3d 


15 


7.354590 


-0.866 


E 


7.39 


2.50 


IPM 


8717.797 






5 


5.932787 






3 


7.354592 


-2.866" 


E 


7.38 






8717.811 






5 


5.932787 






5 


7.354590 


-1.690" 


D 


7.32 






8717.819 






5 


5.932787 






7 


7.354589 


-0.941" 


D+ 


7.39 






8736.012 


3d 


3d 


15 
12 


5.945915 
5.945914 


7/ 


F° 


28 
7 


7.364755 
7.364755 


-0.351 


E 


8.21 


3.00 


IPM 


8736.000 






5 


5.945913 






7 


7.364755 


-3.210<= 










8736.014 






7 
15 


5.945915 
5.945915 




3po 


7 
21 


7.364755 
7.364755 


-1.788" 


E 








8736.000 






5 


5.945913 






5 


7.364755 


-1.786" 


E 








8736.000 






5 


5.945913 






7 


7.364755 


-0.890'' 










8736.014 






7 


5.945915 






5 


7.364755 


-3.337" 


E 








8736.014 






7 


5.945915 






7 


7.364755 


-2.170'' 










8736.014 






7 


5.945915 






9 


7.364755 


-0.725" 


C 








8736.024 






3 


5.945917 






5 


7.364755 


-1.056" 


D+ 









' iHirata fc Horaguchil lll995l') . '^ iKurucz fc Belli lll995l '). 



Table 5. Mgl W/W* values for logg = 1.5. Values outside [0.9, 1.1] range are boldfaced. Blanks mean that W <1 m/ 





[Fe/H] 
%a [K] 




-3 






-2 











A[A] 


3500 


4250 


5250 


3500 


4250 


5250 


3500 


4250 


5250 


4571.095 




0.90 


0.34 


1.46 


0.98 


0.69 


0.27 


0.99 


0.97 


0.73 


4167.271 
4702.990 
5528.403 




0.95 
1.18 
1.37 


0.63 
0.82 

0.96 


0.57 
0.65 
0.71 


1.01 
1.08 
1.19 


0.96 
1.03 
1.11 


0.85 

0.95 
1.04 


0.97 
0.97 
0.97 


0.94 
0.93 
0.92 


0.97 
0.99 
0.97 


8806.751 




1.39 


1.15 


0.91 


1.22 


1.18 


1.24 


0.93 


0.86 


0.95 


5167.320 
5172.683 
5183.603 




0.95 
0.94 
0.94 


0.87 
0.82 
0.79 


0.87 

0.95 
0.97 


0.91 
0.90 
0.90 


0.85 
0.83 
0.82 


0.99 
0.99 
0.97 


0.93 
0.93 
0.93 


0.84 
0.83 
0.83 


0.85 
0.82 
0.81 


4730.028 
5711.086 




1.39 

1.48 


0.79 

0.86 




1.08 
1.22 


0.94 
1.04 


0.77 
0.83 


0.98 
1.01 


0.92 
0.94 


0.89 

0.94 


11828.18 




1.47 


1.23 


0.96 


1.26 


1.24 


1.32 


0.98 


0.87 


0.99 


8923.563 
7657.599t 




1.30 

1.75 


0.85 




1.20 
1.48 


1.28 


0.68 
0.89 


0.96 
1.12 


0.89 

1.05 


0.86 

0.98 


10312.52 
















0.97 


0.95 


1.01 


8997.147 
8473.688 










1.19 


1.09 




1.02 
1.00 


1.00 
0.97 


1.08 
1.03 


(8715.255) f 












0.18 




1.01 


1.10 


1.22 


8736.012 




1.45 


0.82 




1.08 


1.03 


0.56 


1.36 


1.28 


1.23 



t Note that this line is sensitive to the surface gravity. 



Table [S] and Fig. lAl[) . At lower metallicities there is a strong 
variation with TcS and W/W can vary between 0.4 and 1.9, 
which corresponds for weak lines to NLTE abundance cor- 
rections between +0.4 and —0.3 respectively. Their W/W* 
are not sensitive to the surface gravity except for the line at 
8806 A which can have 0.1 dex of difference in abundance 



correction between log g — 1 and log g = 2 at 4500 K and at 
solar metallicity. 



The Mgl AA 4730, 5711 and 11828 lines also come from 
the 3p ^P° level and reach ^ 5 terms. The W/W* of Mgi 4730 
and 5711 A lines have similar trends (0.9 ^ W/W* ^ 1.1 for 
— 1 ^ [Fe/H] ^ +0.5, stronger NLTE variations otherwise) 
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Table 6. Comparison of Mgl abundance corrections with lShinianskava et alj JgOOOl) for five lines as functions of atmospheric parameters. 
We produce NLTE corrections also taking into account the hydrogen collisions by Drawin's formula. The Mg abundances of our results 
follow the enhancement of o-elements of the model atmospheres. Blanks are used when Eq.jGj is not applicable for our NLTE corrections 
or when W is below 1 mA, and when NLTE corrections are not provided bv Shimanskava et al.i (,2000) . 



TW 



Toff [K] 

logs 

[Fe/H] 



4500 
1 



NLTE correction 



A[Mg/H] 

4571.095 A[Mg/H] 0.1 

A[Mg/H]t 0.1 



0.03 -0.01 



0.36 
0.24 
0.14 



5000 
2 



0.08 



0.11 
0.06 
0.08 





A[Mg/H] 





-0.11 


0.09 


;.oi2 


A[Mg/H] 


0.1 


-0.01 


0.00 




A[Mg/H]t 


0.1 


-0.05 


-0.05 



-0.09 

0.00 

-0.05 



0.21 
0.00 



0.48 
0.26 





A[Mg/H] 





0.05 




0.09 






0.13 


4702.990 


A[Mg/H] 


0.1 


-0.01 




0.00 






0.00 




A[Mg/H]t 


0.1 


0.01 


-0.11 


0.02 


0.00 


0.14 


0.14 




A[Mg/H] 





0.20 


0.25 




0.18 






5183.603 


A[Mg/H] 


0.1 


0.10 


0.05 




0.10 


-0.01 






A[Mg/H]t 


0.1 




0.02 


-0.29 




0.00 


-0.14 




A[Mg/H] 







-0.01 


0.04 




0.02 


0.10 


5711.086 


A[Mg/H] 


0.1 




0.00 


0.00 




0.00 


0.00 




A[Mg/H]t 


0.1 


-0.03 


-0.03 




-0.03 


0.12 





t 



Shimanskaya (private communication) 



while the W/W* of the 11828 A line is very similar to the 
W/W* of the 8806 A line. 

The Mgi 8923 A line is not very sensitive to the NLTE 
efTects except for models with Teff > 4500 K and —2 ^ 
[Fe/H] ^ -1 for which we can have A[Mg/H] = +0.15 dex 
of abundance correction for the hottest models. 

The weak Mgt 10312 A line is not sensitive to NLTE 
effects since abundance correction is less than ±0.05 dex. 

The 7657 A line is the counterpart of the 8923 A line in 
the triplet system. This line is in the saturated part of the 
curve of growth at solar metallicity and becomes weak for 
[Fe/H] ^ —1. There is a systematic negative NLTE abun- 
dance correction for models with Tcff ^ 4000 K that can 
reach —0.3 dex. 

The components of the Mgi b triplet follow the same 
NLTE trends. When these components are strong with 
Tefi > 4000 K and with -1 sj [Fe/H] ^^ -H0.5, the NLTE 
abundance corrections are positive and vary between -1-0.1 
and +0.3 dex for the most metal-rich models. For models 
with Tcff ^ 4500 K and with -4 ^ [Fe/H] ^ -3, strong 
variations exist for W/W that can produce NLTE abun- 
dance correction as large as +0.3 dex at [Fe/H] = —4 in 
the linear regime. The sensitivity of the surface gravity to 
the NLTE effect is maximum for the intermediate Teff and 
for the most metal-rich models, and can reach +0.1 dex for 
models between logg = 1 and 2. 

We compare our results with those of 
IShimanskava et al.l l|200d ) for which grids of stellar 
parameters overlap. They provided us NLTE corrections 
for several lines (Shimanskaya, private communication). 
They take into account inelastic collisions with hydrogen 
through the use of Drawin's formula with a scaling factor of 
Sh = 0.1. Therefore, we added hydrogen inelastic collisions 
with the same scaling factor for six model atmospheres 



(Tcff = 4500 and logg = 1) and (Teff = 5000 K and 
logg = 2) with [Fe/H] = 0, -2, and -3. Results for five 
hues (AA 4571, 4702, 5183, 5711 and 8736) are presented 
in Table [6l We formulate three remarks. Firstly, the use 
of Drawin's formula reduces strongly NLTE effects, even 
with a small sca ling factor. Second l y, our NLTE effects 
are smaller than IShimanskava et al.l ([20001) when we use 
Drawin's formula. Thirdly, one should be cautious with 
this comparison because we used different model atoms, 
model atmospheres, atomic data and codes. We noticed 
that we have a strong abundance correction for the in- 
tercombination resonance 4571 A line that can reach 0.5 
dex even if we treat electronic collision of this line not 
using the IPM based on the low log gf value (—5.688), 
but using a collis io nal s trength Qij of 1. Results from 
IShimanskava et al.l ([20001) give lower NLTE abundance 
correction for this line, maybe because they used isotopic 
shifts. The 5711 A line is only affected by NLTE effects 
at lower metallicities. The subordinate 5183 A line of the 
Mgi b triplet is strongly affected by NLTE effects even at 
solar metallicity with (0.1 dex) and without (0.2 dex) the 
use of Drawin's fornmla. 

The Gaia lines 

Five Mgi lines from the triplet system belong to the 
wavelength range of Gaia/RVS. These lines are multiplets 
with poorly known oscillator strengths. The best accuracy 
for a component of the Mgi 8736 A line is C (25%) while 
all the other log gf have an accuracy of D or E, in the 
nist's definition. In order to deliver NLTE effects for these 
multiplets, we use the model atom B with merged levels, 
except for the 8473 A line. To illustrate the fact that the 
SE is not disturbed when introducing the merged levels, we 
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5.946 - 
5.932 - 




3po 



■^D 



3po 



Figure 8. Partial Grotrian diagram for Mgl lines involved in the 
wavelength domain of GaJa/RVS. Mean levels are represented 
by dashed lines with their NIST values. Statistical weights are 
represented on the right side of the fine levels. 



0.2 r 



£ 0.0 r 




4p "P° 
3d ^D 
7d ^D 

7f '-r^ 



ine formation 



Figure 9. Ratio between sum of detailed level populations (n\ + 
^2 +"-3) with average level population (nt) in logarithmic scale 
as a function of standard optical depth. For the 7/ configuration, 
we merge the terms of the 2 multiplicities. We have represented 
domain formations of lines in Gaia/RVS wavelengths for a solar 
MARCS model. 



compare the level populations of the Mgi model atoms A 
and B (so the number of levels drops from 150 at 141 in the 
SE). We compare the sum of populations of fine levels of the 
Mg I model atom A with the population of the merged levels 
of the Mg I model atom B. The logarithm of the formed ratio 
is plotted for each level considered in Fig. (9] The difference 
between populations seems to be important, but when we 
integrate over the depth-scale, the relative error is smaller 
than 0.1% for 4p ^P° and smaller than 0.05% for the three 
higher levels. We conclude that the merging of fine levels for 
4 average levels does not disturb the entire SE of the Mgi 
atom. 

The W/W* for the Mgi Gaia/RVS lines are plotted 
for several Tofi as a function of [Fe/H] for each log g in 
Fig. [BT] The 8473 A line is the least affected by NLTE ef- 
fects but vanishes as soon as the metallicity becomes less 
than —1 dex. The 8715 A Ritz line (representing a fictitious 
fine with the AA 8710, 8712 and 8717 triplet lines) suffers 
from large NLTE effects and it is the most sensitive to sur- 
face gravity. The NLTE abundance correction can vary of 
0.25 dex for T^s = 3500 K, [Fe/H] = -1 between logy = 1 
and 2. The 8736 A line, which is the most visible Mgi line 
in this range, disappears under [Fe/H] = —2. This W/W is 
larger than one for a large range of stellar parameters and 
can reach 1.4 at solar metallicity. 

In our knowledge, there is only N LTE results for the 
Mgi 8736 A line in the literature by IShimanskava et al.l 
l|200Cl) . A comparison for six model atmospheres are pre- 
sented in Table [S] For models at solar metallicity, we are 
in 0.05 dex agreement. But for the model at Tes = 4500, 
logg = 1 and [Fe/H]= —2, we find a positive NLTE abun- 
dance correction (for Sh = 0). This positive NLTE abun- 
dance correction increases with effective temperature and 
surface gravity whereas this NLTE abundance correction is 
negative and constant at —0.05 dex for Shimanskaya (private 
communication). This can be due to the difference in the Mg 
abundance chosen at this metallicity ([Mg/Fe] — +0.4 for us 
versus +0.00 for them) and to the value of Sh- 



4.4 The Cai and Can lines 

4-4- i The Cai lines 

We show the departure coefficients bi for the Cai model 
atom in Fig. 1101 with the same model atmosphere used for 
the Mgi model atom in Figl?! i.e. T^a = 4500, logg = 2 
and [Fe/H] = —3. The ground stage of Cai is also depopu- 
lated by overionization but in a lesser extent than for Mg i. 
In general, the Cai model is less affected by NLTE effects 
than the Mg i model for the atmospheric parameters used in 
this paper. Two mechanisms appear: on the one hand, over- 
photoionization depopulates lower levels and on the other 
hand, coUisional recombination and photon suction overpop- 
ulate levels close to the ionization stage. 

The Cai 4226 A resonance line shows that its W/W* 
is strongly dependent of the three atmospheric parameters 
(see Table [3 and Fig. \M^. For aU the models, the NLTE 
abundance corrections are positive with A[Ca/H] ^ 0.03. 
This line is very sensitive to the surface gravity except for 
Tcff > 5000 K and [Fe/H] sj -2. The NLTE abundance 
correction can reach 4-0.2 dex between logg — 1 and 2 for 
the most metal-rich models and -1-0.5 dex for the coolest and 
most metal-poor models. 

The Ca i 6572 A intercombination and resonance line is 
also sensitive to the surface gravity but in a lesser extent. 
This line is weaker than the 4226 A one and shows that 
W/W* depends on the logp only for models with Tcff ^ 
4000 K. 

The Cai AA 5512, 5867 and 7326 lines share the same 
lower level 4p ^P° . For these lines, 0.5 ^ W/W* ^ 1.2 except 
for the coolest and most metal-poor model. We can adopt 
that they are formed under LTE condition in the range of 
— 1.5 ^ [Fe/H] ^ 0. NLTE abundance correction can be as 
large as +0.3 dex for [Fe/H] ^ -1.5. For Tos < 4000 K and 
[Fe/H] < —3, the 5512 and 7326 A lines become sensitive to 
the surface gravity. 

For the Cai red triplet AA 6102, 6122 and 6162 sub- 
ordinate lines, the NLTE abundance correction are positive 
for almost all the models. Large deviations from LTE ex- 



A grid of NLTE corrections for Mg and Ca for evolved stars 15 

Table 7. Cal W/W* values for logg = 1.5. Values outside [0.9, 1.1] range are boldfaced. Blanks mean that W <1 mA. 
A[A] fFc/Hl -3 -2 



[Fc/H] 
Toff [K] 



3500 4250 5250 



3500 4250 5250 



t Note that this line is sensitive to the surface gravity. 



3500 4250 5250 



4226.727T 


0.54 


0.59 


0.93 


0.64 


0.66 


0.86 


0.84 


0.58 


0.81 


6572.777T 


0.58 


0.35 




0.85 


0.67 


0.51 


0.93 


0.81 


0.75 


5512.980 


1.02 


0.61 


0.52 


1.04 


0.92 


0.76 


0.95 


0.97 


1.04 


5867.559 


1.28 






1.04 


0.92 


0.76 


0.99 


0.87 


0.97 


7326.142 


1.13 


0.71 


0.51 


1.10 


1.08 


0.85 


0.94 


0.97 


1.08 


4425.435 


0.83 


0.59 


0.52 


0.85 


0.90 


0.82 


0.81 


0.85 


0.96 


6102.719 


0.96 


0.64 


0.56 


0.88 


0.99 


0.86 


0.77 


0.83 


0.96 


6122.213 


0.96 


0.77 


0.60 


0.86 


1.00 


0.97 


0.75 


0.80 


0.98 


6162.170 


0.92 


0.82 


0.66 


0.84 


0.97 


1.03 


0.74 


0.76 


0.98 



6161.294 


0.57 


0.40 




0.76 


0.60 


0.58 


0.87 


0.77 


0.83 


6166.438 


0.55 


0.40 




0.73 


0.60 


0.57 


0.87 


0.78 


0.85 


6169.037 


0.55 


0.39 


0.39 


0.79 


0.66 


0.57 


0.86 


0.80 


0.88 


6169.562 


0.59 


0.42 


0.46 


0.83 


0.75 


0.65 


0.84 


0.83 


0.90 



4578.549 


0.49 


0.36 


0.46 


0.77 


0.65 


0.59 


0.88 


0.87 


0.91 


6455.592 


0.56 


0.40 




0.77 


0.62 


0.62 


0.92 


0.82 


0.91 


8525.707 
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0.62 


8583.318 
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Figure 10. Cal departure coefficient as a function of optical 
depth in the continuum at 5000 A for a metal-poor giant stellar 
model (Teff = 4500 K, logg = 2.0, [Fe/H] = -3.0 and [a/Fe] = 
+0.4). We separate levels as functions of their energy Ei in respect 
to the ground state as shown in the legend. 
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Figure 11. Global Cal departure coefficients for different model 
atmospheres of T(,ff = 4500 K, logg = 1 (thin lines), logg = 2 
(thick lines), and [Fe/H] = (solid line), [Fe/H] = -3 (dotted 
line); and of T^s = 5250 K, logg = 1 and [Fe/H]= -1 (dashed 
line). The fourth parameter on the plot is the o-enhancement of 
the model atmosphere. 



ist for Tcff < 4500 K with -1 ^ [Fe/H] sC +0.5, and for 
Tcff > 4500 K with -4 ^ [Fe/H] < -2.5. They are not 
sensitive to the surface gravity. 

The Cai lines at AA 6161, 6166, 6169.0 and 6169.6 be- 
long to the same multiplet. Their W/W behaviours are 
very similar and always less than one, that produces posi- 
tive NLTE correction. The NLTE abundance corrections can 
reach +0.4 dex for [Fe/H] ^ —3. Theses lines are sensitive 
to the surface gravity for models with Toff < 4000 K and 
[Fe/H] < -3. 

The line at 4578 A is also from the same lower term 
than the previous multiplet. The W/W* is very similar at 
low metallicities. We can note that this line is formed in 



LTE for models with -1 ^ [Fe/H] < (W/W*between 0.9 
and 1). 

The Cai 6455 A intercombination line also shares the 
same lower term than the previous multiplet and the 4578 A 
line. The W/W is also similar and always lower than one. 
We notice that the lower deviation from LTE comes near to 
[Fc/H] = -1. 

The Cai Gaia/RYS lines at 8525, 8583 and 8633 A 
are strongly affected by NLTE effects at solar metallicity 
(see Table (Tjl. These lines are weak whatever the stellar pa- 
rameters (NLTE EW ~ 50, 70 and 90 mA at maximum 
respectively). The stronger NLTE effects appear for 8525 
and 8583 A lines and can reach A[Ca/H] — +0.4 dex. For 
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Tcff > 4250 K (respectively T^ff < 4250 K), NLTE effects 
increase (respectively decrease) with decreasing metallicity 
as shown in Fig. IBll 

W e compare our NLTE results with those of iDrakd 
l|l99H ) and found opposite effects concerning the variation 
of overion ization w i th sur face gravity. As emphasized in the 
review of lAsplundl (|2005l ). Drake unexpectedly found that 
the over ionization of Cai decreases with surface gravity 
whic h is expla i ned b y an opacity effect that did not con- 
vince lAsglliildl (|2005l ). For giants stars, we show in Fig. 1111 
th e oppos i te effe cts of overionization, comparing with fig. 6a 
by iDrakd (|l991h . We plot overall Cai departure coefficient 
for model atmospheres of Tofi = 4500 K, logg — 1 and 2, 
and [Fe/H] = and —3. We see that for a given metallicity, 
the overall Cai departure coefficient tends to be closer to 
one with the increase of the surface gravity. This is the op- 
posite of the two models in fig. 6a of Drake (Toff = 4500 K, 
log*; = 1 and 0, and [Fe/H] = 0) where the overionization 
of Ca I increases when decreasing the surface gravity . 

Comparison with the work of iMashonkina et al.l (|2007l ) 
is done even if the grids of stellar parameters do not overlap. 
They provided NLTE abundance corrections for a grid of 
dwarfs and subgiants. We compare results for the Sun and 
for a spherical model at T^g = 5000, logg — 3, [Fe/H] = —2 
and [a/Fe] — +0.4. In order to be the most consistent for 
the comparison, we use our Cai/ii model atom described 
in Sect. 12.3.31 Results are shown in Table [§1 We found a 
better agreement with these authors for the metal-poor star 
rather than for the Sun. Our NLTE abundance corrections 
for the solar Ca i lines are mainly lower and in opposite sign 
compared to them. For the metal-poor star, agreement is 
very satisfactory except for the 4425 A line for which we 
have a halved abundance correction. These agreements are 
found using different geometries f or the model atmosphere s 
(plane-parallel Kurucz models for IMashonkina et al.l (|2007l ) 
against spherical marcs models for us) and different NLTE 
codes, but give confidence for the NLTE computations done 
for the grid of late-type giant and supergiant stars. 



4.4-2 The Can lines 

A part of the results for Ca ii lines are shown in Table [S] 
and in Fig. IA3I For the strong resonance H & K lines 
at AA 3968 and 3933 respectively, W/W is less but close 
than one except for temperatures below 3900 K. Hence, the 
NLTE abundance correction does not exceed 0.1 dex be- 
tween [Fe/H] = -1-0.5 and —3. It can vary between —0.1 
and 0.2 dex at the lowest metallicity of —4 dex. We no- 
ticed a light sensitivity to the surface gravity for all effective 
temperatures. For [Fe/H] ^ —2.5 that can reach 0.1 dex of 
amplitude between logg = 1 and logg = 2 at 3500 K (see 
Tableland Fig. lX3| . 

The lines at AA 8248 and 9854 come from the same 
5p ^P° term but from different fine levels and have W/W* > 
1.1. The deviation from LTE increases with increasing effec- 
tive temperature and can reach —0.1 dex of abundance cor- 
rection. These weak lines are sensitive to the change of the 
surface gravity for T^h ^ 4500 K. These variations can reach 
0.05 dex of amplitude between logg = 0.5 and logg = 2. 

The NLTE abundance correction for the 8912 A line can 
vary between —0.15 dex at solar and subsolar metallicities 
and -0.3 dex at [Fe/H] = -3.5 and T^a = 4500 K. 




Figure 12. EW ratio of Ca II IR 8498 A line as a function of effec- 
tive temperature T^g and metallicity [Fe/H], for a given surface 
gravity logg = 1.5. 
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Figure 13. NLTE/LTE EW ratio of Can IR 8662 A line as 
function of metallicity [Fe/H] and logg, for a given effective tem- 
perature T^ff = 5000 K. Each colour of the W/W represents a 
surface gravity (red, green, blue and cyan for logg = 0.5, 1.0, 1.5 
and 2.0 dex respectively). Dots represent computed ratios while 
full lines represents our polynomial fit of the NLTE/LTE EW ra- 
tios. The dashed curve stands for the W/W * {8662 A) polynomial 
fit computed bv lStarkenburg et al.l l|201Cl ). 



The Can IR 11949 A line shows an NLTE sensitivity 
to the surface gravity (as seen in Fig. IA3|l that increases 
with the effective temperature. W/W* > 1.2 for all of the 
atmospheric parameters. 

The NLTE/LTE EW of each component of the CaT 
follows the same trend. An example of the W/W* trend of 
the CaT is shown for the 8498 A line, which has the most 
pronounced NLTE effects among the triplet, in Fig. [12] for a 
surface gravity of logg = 1.5. NLTE effects on the CaT are 
weak at solar and sub-solar metallicities for all the effective 
temperatures. For lower metallicities, strong deviations from 
LTE appear. As shown in Fig. [12 [13] and IH the W/W* 
trends of the CaT as a function of the metallicity can be 
separated in two regimes: 
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Table 8. Call W/W* values for logg = 1.5. Values outside [0.9, 1.1] range are boldfaced. Blanks mean that W <1 mA. 





[Fc/H] 
Teff [K] 




-3 






-2 











A[A] 


3500 


4250 


5250 


3500 


4250 


5250 


3500 


4250 


5250 


3933.663 
3968.468 






0.97 
0.96 


0.95 
0.95 


0.98 
0.97 


0.96 
0.97 


0.98 
0.98 


0.99 
0.99 


0.99 
0.99 


0.99 
0.99 


0.99 
0.99 


8498.018* 
8542.086* 
8662.135* 






1.13 
1.13 
1.11 


1.06 
1.02 
1.02 


1.23 
1.11 
1.14 


1.02 
1.01 
1.00 


1.00 
1.00 
1.00 


1.07 
1.04 
1.05 


0.95 
0.95 
0.95 


0.94 
0.95 
0.95 


0.98 
0.98 
0.98 


8248.791* 
9854.752* 










1.12 




1.17 


1.24 
1.24 


1.13 
1.12 


1.23 
1.22 


1.26 
1.27 


8912.059 








1.91 


1.54 




1.55 


1.64 


1.27 


1.37 


1.38 



11949.74* 1.69 1.57 1.26 1.51 1.80 1.27 1.43 1.51 

t Note that this line is sensitive to the surface gravity. * Note that this line is sensitive to the surface gravity for the hottest models. 



Table 9. Comparison of Cal abundance corrections with iMashonkina et al.l 1 12007^ for the Sun and for a metal poor sub-giant star 
without hydrogenic collisions (Sh = 0)-W represents our NLTE computed EW and W/W* represents our NLTE/LTE EW ratios. When 
Eq. (O is not applicable, blanks are used in A[Ca/H] column. Blank in A[Ca/H]* column denotes a NLTE EW less than 5 mA. 



Sun 

A [A] W [mA] W/W* A[Ca/H] A[Ca/H]* 



Tefi = 5000, logg = 3, [Fe/H] = -2 
W [mA] W/W* A[Ca/H] A[Ca/H]T 



Cal 



4226.727 


1318 


0.988 


0.01 


0.07 


543 


0.766 


0.23 


0.24 


4425.435 


174 


1.005 




0.04 


84 


0.830 


0.08 


0.17 


5512.980 


96 


1.056 


-0.01 


0.03 


20 


0.766 


0.12 


0.10 


5867.559 


21 


1.018 


-0.01 


0.06 


2 


0.765 


0.12 


0.08 


6162.170 


260 


1.037 




0.01 


128 


0.990 




0.00 


6166.438 


65 


1.013 


-0.01 


0.07 


9 


0.557 


0.25 


0.30 



Can 



8248.7907 
8498.0180 



68 
1244 



1.184 
1.004 



-0.15 
0.00 



-0.11 
-0.02 



5 
731 



1.153 

1.087 



-0.06 
-0.07 



-0.08 



* IMashonkina et al.l l l2007f) 
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Figure 14. Departure coefficients for fine levels implied in the 
CaT lines for comparison between a metal rich model atmosphere 
(bold curves) and an extremely metal poor model (thin curves). 
The two model atmospheres have Tgff = 5000 K and log g = 2, but 
different global metallicities ([Fe/H]= -1-0.25 and —3.5) and dif- 
ferent a-enhancement factors ([a/Fe]= and +0.4 respectively). 



• the W/W* is a linear function of the metallicity [Fe/H] 
when W/W* s; 1; 

• the W/W* increases strongly with decreasing metallic- 
ity when W/W* > 1 and can reach 1.45. 

Roughly speaking, the CaT lines are formed in LTE 
{W/W* = 1 ± 0.05 which corresponds to A[Ca/H] = 
0.00 ± 0.04 dex) if -2 sj [Fe/H] sj for Tch and logg in 
our grid. The CaT lines cannot be considered in LTE when 
[Fe/H] < -2, and when [Fe/H] > for T^s < 4000 K. The 
NLTE effects decrease with effective temperature and pass 
through a minimum at ~ 4000 K to increase again at lower 
Tcff , as shown in Fig. 1121 This effect may be due to the in- 
crease of molecular continuous opacities in marcs models 
when the effective temperature decreases. 

The behaviour of the W/W* with varying [Fe/H] of the 
CaT can be explained with the variations of the departure 
coefficients presented in Fig. [TJ] As the stimulated emission 
can be neglected for optical and near IR lines, the line source 
function relative to the Planck function Sl/B,y follows the 
departure coefficient ratio bj/bi (i and j are the lower and 
upper levels respectively). The behaviour of the W/W* de- 
pends on the values of bi 
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Table 10. C o mpar ison of our W/Wj with those of 
|j0rgensen et alj 1 119921 ) for the same stellar parameters. The 
NLTE/LTE EW ratio is defined as W/W} = VK(8542A + 
8662A)/W'*(8542A- 



Toff [K] 


log 3 
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W/W* 
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0.95 




1 
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0.95 


0.94 
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-1.00 
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0.99 






0.00 


0.98 


0.99 




1 
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0.99 


0.99 
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-1.00 


0.99 


0.98 






0.00 


0.98 


0.99 




2 


-0.50 


0.99 


0.99 






-1.00 


1.00 


0.99 












|j0reensen et alJ 


Jl992l) 





tion between them. In the metal-rich model (in Fig. ll4p . the 
levels are over-populated due to overionization of Cai and 
bi < bj that implies Sl > Bv, and the emergent intensity 
is strengthened in the line compared to LTE intensity. Thus 
the EW in NLTE is lower than the EW in LTE. This ex- 
plains why W/W* ^ 1 for models with -1 <[Fe/H]< -1-0.5 
in Fig. 1131 For the metal-poor model in Fig. 1141 the mech- 
anism is the opposite due to the change in the relative val- 
ues of bi and 6-,. For metallicities lower than — 1, fe^ > bj if 
Tsooo ^ —1.7 and then Sj, < B^. This implies a reduction 
of the emergent intensity in the line and then a larger value 
of EW in NLTE relative to the LTE. This explains why 
W/W* > 1 for models with [Fe/H]^ -1 in Fig. [H The 
large effects at lower metallicity is due to a larger extent of 
log \bj /bi\ compared to solar and sub-solar metallicities. 

The fine structure of the levels follows the same trends 
but level population differences increase with decreasing 
logrsooo- The Can IR 8498 A line has the greater NLTE ef- 
fects relative to the two other lines because it has the largest 
amplitude on the deviation between the lower coefBcient bi 
(3d ^D^/2, blue dotted line) and the upper coefBcient bj 
(4p ^P3°/2, red dashed line) as seen in Fig. [14] The CaT lines 
are domi nated by their wings even at [Fe/H] = —3 as shown 
in fig. 1 of lStarkenburg et al.l (|2010l ). These wings are formed 
in the deep photosphere in LTE conditions. Therefore the 
EW are weakly influenced by NLTE effects at solar metal- 
licity and moderate for most metal-poor stars for which the 
correction can become very large (20-30% at [Fe/H]= —3). 
These results are in quite good agreement with those of pre- 
vious investigations with the same stellar parameters (e.g. 
IJoreensen et al.lll99l lAndretta eralll2005h . 

We note that even if NLTE line formation has almost 
no influence on Ca abundance determination using the CaT 
at solar or sub-solar metallicity, it has an impact on the 
detection of stellar activity. Indeed, the central depression 
of the CaT is often used as an indic ator of activity, (e.g. 
lAndretta et al.ll2005l . rBusa et al.l2007h . The detail modelling 

of the line core formed in NLTE is therefor e crucial. 

Comparing our results with those of |j0rgensen et al.l 

IJ1992I ). we find that they are in very good agreement. They 
combine the two strongest lines of the CaT, noted here 
W/W J (8542 and 8662 A) and compute NLTE effects for 



a large set of gravities (from logg = 0.0 to 4.0) but for 
a small metallicity range (from —1.0 to 0.2). We give the 
comparison between our result for the combined lines ratio 
W/W} in Table [TO] Deviation in NLTE effect with these 
authors are less than 2%. The NLTE EW increases with in- 
creasing metallicity (— 1 < [Fe/II]< -1-0.25) but our EW are 
larger by a factor of 2 with those of j0rgensen et al. ones 
at [Fe/H] = —1. This may come from the fact that their 
model atmospheres are not a-enhanced at lower metallic- 
ities. Note also that they do not have the photoionization 
cross-sections from the TopBase and ABO theory. Other au- 
thors computed NLTE effects for the CaT but for late-type 
dwarf stars (e.g. lAndretta et al.ll2005r ). They found values of 
W/W* always larger than 1 reaching 1.35 at lower metallic- 
ity ([Fe/H]= -2). 

We also compare d the N LTE effects on CaT computed 
bv lStarkenburg et al.l l|201Cl ) and found a clear discrepancy 
at low metallicity. They provided a 2D polynomial fit of 
W*/W = /(Toff, [Fe/H]) for the lines at 8542 A and 8662 A, 
considering that these ratios are insensitive to the surface 
gravity (1 ^ logp ^ 2). We plotted in Fig. [T3l and IMl their 
polynomial fit for the 8662 A line. The discrepancies be- 
tween our W/W* and theirs, at a given Tcff, increase with 
decreasing metallicity and with decreasing surface gravity. 
At [Fe/H] = —4, the deviation is about 10% for log g = 2 and 
about 18% for loggr = 1. These differences may come from 
the fact that they used a different geometry (plane-parallel) 
in their study. We checked that the use of the scaling fac- 
tor Sh ~ 1 for the collisions with neutral hydrogen cannot 
compensate this difference. Note that they not clearly in- 
dicate the value of Ca abundance that they used for the 
computations of the W/W* . 



Since Starkenburg e t al.l (|2010|) used the model atom of 

[Mashpnkina ot al. (2007|), we investigate the effects of the 
different atomic physics injected in the two model atoms of 
Ca i/ii. We consider a stellar model (Toff — 5000 K, log g — 3, 
[Fe/H] = —2), which is outside of our grid, but used in 
iMashonkina et al.l ((20071) to compare with. Our predictions 
for the two Can lines (8248 and 8498 A) in this particular 
case are shown in Table [S] We note a difference of 0.04 dex 
for the weak line at 8248 A for the Sun and a very good 
agreement (better than 0.01 dex) for the Can IR 8498 A 
line for the meta l-poor model atmo s ph ere. Therefore, the 
discrepancy with IStarkenburg et al.l (|201Q ') does not come 
from the model atoms. It must be investigated in more de- 
tails since the calibration between the CaT and [Fe/H], as 
defined in their paper, may be affected for RGB metal-poor 
stars. 

Interesting opposite NLTE trends appear between Cai 
and Can lines. In general, W/W* for Cai lines are smaller 
than 1.0, whereas W/W* for Can lines are larger than 
1.0, whatever the stellar parameters. We note that the 
NLTE effects are anti-correlated between the Cai AA 6102, 
6122, 6162 triplet lines and the CaT lines in metal-poor 
model atmospheres. This may have a large impact on the 
Ca abundance analysis in metal-poor red giants; such re- 
mark has already been mentio ned in the Gaia context 
JRecio-Blanco fc Theveninll2005l ). 
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Table 11. Fit coefficients of NLTE/LTE EW ratio [W/W*]flt for lines in the GaJa/RVS wavelength range. The fits are valid for 
3500 i£ Teff ^ 5500 K, 0.5 ^ logg ^ 2 and [Fe/H]n,in ^ [Fe/H] ^ +0.5. The Mg and Ca abundances follow the a-enhancement given 
in Eq. l(7)l. The blanks denote coefficients related to discarded terms, emax and emcan are given (in %) as the maximum and the mean 
absolute difference between [Vy/Vy*]flt and W/W*. 
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4.5 Uncertainty estimations 

Uncertainties of W/W* rely on the atomic parameters and 
on the assumptions of the stellar atmosphere modelling 
(viz ID, static, NLTE radiative plane parallel transfer with 
MULTi code using ID, static, LTE, spherical marcs model 
atmospheres). We computed W/W* varying atomic param- 
eters in order to estimate uncertainties of the input atomic 
physics. We selected four model atmospheres (Tch = 5000 K, 
log 5 = 2, [Fe/H] = 0, -2; and T^a = 4500 K, log g = 1, 
[Fe/H] = 0, —2 and three lines (one resonance, one subordi- 
nate, and one allowed lines) for each model atom (AA 4571, 
5183, 8806 for Mgi, AA 6162, 6166, 6572 for Cai, and 
AA 3933, 8248, 8498 for Can). The Mgi 4571 A and the 
Cai 6572 A lines are the equivalent intercombination res- 
onance lines for Mgi and Cai model atoms. For each line, 
the uncertainties of three atomic parameters were tested: the 
oscillator strengths, the photoionization cross-sections and 
the effective coUisional strengths with electrons. We empha- 
size that we performed these uncertainty estimations over a 
hundred of runs (3 among 45 lines over 4 among 453 model 
atmospheres) to give a rough idea of the impact of the un- 
certainties on atomic data. 

We tested for each model atoms the change of 
the oscillator strengths of the three lines by ±50%. 
The changes on the W/W* for the three lines are 
less than 8%. For instance, with the model at T^s = 
5000 K, log g = 2, [Fe/H] = -2, we obtained 



W/W* {Mgi 8806 A) = 1.35J:o.o6. W^/W*(Cai 6572 A) = 
0.53to;o2 and W/W*(Caii 8498 A) = im^o^l- We remind 
that the accuracy of the log gf used and given in Table [3] 
and 13] varies between 2% and 50%. Therefore, these uncer- 
tainties can be seen as an upper limit. We noticed that the 
W/W* are more affected by changes of oscillator strength 
values for the metal-poor model atmospheres. 

Varying the amplitudes of the photoionization cross- 
sections does not affect the lines of model atoms in the same 
way. We changed by a factor ±2 the amplitude of the Top- 
Base photoionization cross-sections of the lower levels of the 
lines considered. The Mgi model atom is the most affected 
with variations on the W/W* that can reach -1-8% (for in- 
stance W/W* {Mgi 5183 A) = 0.79lo;o3 for the model with 
Tcff = 4500 K, log 5 = 1, [Fe/H] = 0). Calcium is less 
affected, with 1% of variations of W/W* for the Cai and 
less than 1% for t h e Ca n model atoms. As emphasized by 
iMashonkina et al.l ((20071), the changes of the photoioniza- 
tion cross-sections mainly affect the minority species (Cai 
in the range of stellar parameters considered here). How- 
ever, we found that the absolute corrections on the W/W* 
increase with decreasing metallicity, due to variations on the 
amplitude of photoionization cross-sections. 

Finally, we tested the variations of the amplitude of 
the collisional strength with electrons. Using the IPM ap- 
proximation, the change by a factor of ±2 of the effective 
collisional strength slightly affects the W/W* results (less 
than < 2%). 
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4.6 Polynomial fit of the W/W* for the Gaia lines 

Due to the importance of the CaT and the Mg i 8736 A for 
the Gaia/RVS mission, we provide a fit of the W/W pre- 
sented in the previous section as a multivariable polynomial 
for the five Mgi lines, the three Cai lines and the CaT lines. 
The expression of the fit is up to the third order in stellar 
parameters: 



(8) 
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with: 

X = (Tcff - 4375) /875 
y = (logg-1.25)/0.75 



(9) 



([Fe/H] - [Fe/H]c)/(0.5A[Fe/H]) 



the reduced and centred variables that lead to a variation in 
the [—1, 1] range for our ranges of atmospheric parameters. 
A [Fe/H] represents the metallicity range and [Fe/H]c is the 
median metallicity, which depends on the line considered. 
As [Fe/H]max = +0.50 in all cases, the median metallicity 
and the metallicity range are expressed by: 



[Fe/H]c = -[Fe/H]^in + 0.25 



A [Fe/H] = [Fe/H]„ 



[Fe/H]n,in 



where [Fe/H]min is specified in the Table 1111 As seen in 
Figs. |BT] and IHl the NLTE/LTE EW ratios have different 
behaviours in stellar parameters for each li ne considered. A 
second ord er formula (JAndretta et aLlbOOSl *) or modified sec- 
ond order (jStarkenburg et al.ll201Gl ) is not enough to fit the 
dependence, especially for [Fe/H] at low metallicity (^ —3). 
In order to find the coefficien ts am {m = 0, 19), we use the 
LSQ package of lMillen (|l992l ) and we rejected W/W* values 
for which W < 1 mA. We present the results in Table [11] 
for the Mg i, Ca i and Ca ii Gaia/RVS lines. We noticed that 
some terms may be discarded without modifying the accu- 
racy of the fit: the related coefficients are represented by 
blanks in Table 1111 For instance, we see that for Mg i and 
Cai lines, the fits are insensitive to the xy^ and y^ terms; 
for the CaT, the fits are insensitive to the xy, y^ and xyz 
terms. 

Examples of the fits are shown in Fig. [13] for the Ca ii 
IR 8662 A and in Figs. [Bl] and [B2] for the other Gaia/RYS 
lines. For the Gaia/RVS Mgi and Cai lines we restrict the 
range of metallicities as indicated in Table [TT] since below 
the lower metallicities the line is too weak (W ^ 1 mA) 
to be considered. These fits can be used to estimate the 
NLTE/LTE EW ratios with an accuracy better than 10%, 
in the range of the stellar parameters considered. The largest 
deviations of the fits appear for the most metal-poor model 
atmospheres for the Mg i 8736 A and the CaT lines. Hence, 
when requested accuracy is greater than emcan (given in % in 
Table [11}, our fits can be used with confidence. Otherwise, 
it is advised to directly use W/W in the electronic tables. 



5 CONCLUSION 

We have performed NLTE computations for the Mgi, Cai 
and Can model atoms in late-type giant and supergiant 
stars. We provide NLTE/LTE EW ratios for a grid of 
453 MARCS model atmospheres (3500 ^ Teg ^ 5250 K, 
0.5 s; \ogg 4: 2.0 dex and -4.00 ^ [Fe/H] sC -1-0.50 dex) 
in electronic forma^^l. The model atoms are based on the 
assumption that we do not take into account inelastic col- 
lisions with neutral hydrogen since realistic quantum me- 
chanical collisional cross-sections are still unavailable. We 
used a formulation for electronic collisions that underesti- 
mate the collisional rate. Such underestimation induced an 
upper limit on the NTLE/LTE EW ratios especially for the 
Mgi lines in the Gaia/RVS wavelength range. The use of fine 
structure in the model atoms do not affect strongly NLTE 
results since the W/W of the components of a multiplet 
are very similar but permit a consistent representation of 
the physics. This work will be extended to late-type main 
sequence stars when collisional cross-sections with hydrogen 
are available. 

The main conclusions for the lines outside the 
GaJa/RVS wavelength range are as follows. For the Mgi 
and Cai lines, the assumption of LTE underestimates the 
Mg and the Ca abundances and gives mainly positive NLTE 
abundance correction. For the Can lines, the assumption of 
LTE overestimates the Ca abundance and gives mainly neg- 
ative NLTE abundance corrections. However, most of the 
Mgi lines show W/W > 1 that can reach 2 for the coolest 
models with an increase of the sensitivity to the surface grav- 
ities. The Mgi b and Cai red triplets show W/W* < 1 that 
can reach 0.5 for the most metal-poor and the hottest model 
atmospheres. The NLTE effects for the Mgi 4571 A and the 
Cai 6572 A intercombination and resonance lines are very 
sensitive to the metallicity and to the surface gravity. The 
NLTE effects on Can H&K lines are negligible except at 
[Fe/H] = -4. 

For the Gaia/RVS lines, NLTE computations give the 
following trends. The very weak 8473 A line is mainly formed 
in LTE whereas the Mgi 8736 A and the Mgi AA 8710, 8712 
and 8717 triplet lines are mainly formed in NLTE with a 
strong sensitivity to the surface gravity for the triplet. The 
weak Ca i triplet lines are mainly formed in NLTE but vanish 
as soon as [Fe/H] < 1. The famous CaT lines are mainly 
formed in LTE if [Fe/H] ^ -2. The NLTE effects increase 
with a decrease of the metallicity and with an increase of the 
surface gravity. We show that the H^/W*(CaT) can increase 
by 20% for log g varying from 0.5 to 2.0. For convenience, we 
provide a polynomial fit computed for the Gaia/RVS lines: 
the Mgi AA 8473, (8715) triplet, 8736 muhiplet, the Cai AA 
8525, 8583 and 8633, and the CaT lines. The polynomial 
can be extensively used by the automatic tools dedicated to 
the analysis and extract chemical abundances for millions 
of stars in the context of the large surveys as, for example, 
Gaia and RAVE. 



^^ available on the Vizier/CDS database 
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APPENDIX A: W/W* FOR THE LINES 
OUTSIDE OF THE GAM/RVS WAVELENGTH 
RANGE 

We present the results of the NLTE computations as the evo- 
lution of the NLTE/LTE EW ratios W/W* as a function of 
the metallicity [Fe/H]. Each row represents a selected line. 
Each panel of a row represents the W/W* for a given ef- 
fective temperature (3500, 4250 and 5250 K). Each colour 
of the W/W* in a panel represents a surface gravity (red, 
green, blue and cyan for logg = 0.5, 1.0, 1.5 and 2.0 dex re- 
spectively). The variations of W/W* are not plotted when 
W < 1 mA. We notice that marcs model atmospheres are 



missing for logg = 0.5, T^« G [3500, 4500], and [Fe/H] = -4; 
and for \ogg = 0.5, T^a £ [3500, 3900], and [Fe/H] = -3.00. 
The dotted line stands for no deviation from LTE. The 
dashed lines stand for a variation of ±5% of the W/W* 
(corresponding to a variation of ±0.04 dex of the A[E1/H]) 
in case of strong regime of the curve of growth theory. The 
solid lines stand for a variation of ±10% of the W/W* (cor- 
responding to a variation of ±0.04 dex of the A[E1/H]) in 
case of weak regime of the curve of growth theory. 

We show the most representative Mg i lines in Fig. lAll 

• the 5183 A W/W* represents the Mgi b triplet lines at 
AA 5167, 5172 and 5183; 

• the 5711 A W/W* represents the similar behaviours of 
the A A 4730, 5711 and 11828 hues; 

• the 7657 A W/W*; 

• the 8806 A W/W* represents the similar behaviours of 
the AA 4167, 4702, 5528 and 8806 lines that are issued from 
the 3p ^P°; 

• the 8923 A W/W* that is the counterpart of the 7657 A 
in the singlet system; 

• the 10312 A W/W* that is essentially formed in LTE. 

We show the most representative Ca i lines in Fig. IA2I 

• the 4226 A W/W* represents the NLTE behaviour of a 
resonance line; 

• the4578 A W/W*; 

• the 5512 A W/W* represents the similar behaviour of 
the AA 5512, 5867 and 7326 lines; 

• the 6122 A W/W* represents the similar behaviour of 
the Cai red triplet lines at AA 6102, 6122 and 6162; 

• the 6166 A W/W* represents the similar behaviour of 
the multiplet lines at AA 6161, 6166, 6169.0 and 6169.6; 

• the 6572 A W/W* represents the NLTE behaviour of 
an intercombination resonance line. 



We show all the selected lines for the Can lines in Fig. 



APPENDIX B: W/W* FOR THE LINES IN THE 
G^M/RVS WAVELENGTH RANGE 

We present the results of the NLTE computations and the 
polynomial fits as the evolution of the NLTE/LTE EW ra- 
tios W/W* as a function of the metallicity [Fe/H] for the 
Mg I, Ca I and Ca ii lines in the Gaia/RVS wavelength range 
[8470, 8740] A in Fig. [Bl] and [H Each row represents an 
observed line, except for the quoted row at 8715 A standing 
for the Ritz wavelength of the Mgi triplet at AA 8710, 8712 
and 8717. Each panel of a row represents the W/W* for a 
given effective temperature (3500, 4250 and 5250 K), except 
for the Can 8542 and 8662 A lines for which we used the 
following temperatur es (3500, 4250 and 5 000 K ) in order to 
add polynomial fit of IStarkenburg et al.l l|2010l ) on the pan- 
els. Each colour of the dots and the polynomial fits W/W* 
in a panel represents a surface gravity (red, green, blue and 
cyan for logg = 0.5, 1.0, 1.5 and 2.0 dex respectively). 

This paper has been typeset from a TJjX/ ETJjX file prepared 
by the author. 
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Figure Al. W IW* for the selected Mgl lines as function of the stellar parameters (see Appendix A for details). 
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Figure A2. W/W* for the selected Cal lines as function of the stellar parameters (see Appendix A for details). 
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Figure A3. W/W for the selected Call lines as function of the stellar parameters (see Appendix A for details). 



26 T. Merle, F. Thevenin, B. Pichon, L. Bigot 



1.0 
0.8 
0.6 



^ - 

S 1.0 :■■ 

I 0.8^ 

> - 

•^ 0.6 : 

> 0.-1 - 




0.5 0.0 0.5 1.0 

[Fc.'H] 



Figure Bl. W/W* for the Gaia/RVS Mgl and Cal lines as function of the stellar parameters (see Appendix B for details). 
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Figure B2. W/W* for the Gaia/RVS Call lines as function of the stellar parameters (see Appendix B for details). 



